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1.0  INTRODUCTION 


This  document  contains  a  bibliography  of  reports,  papers  and  books 
that  deal  directly  with  Prony's  method.  No  attempt  has  beri  made  to 
Include  papers  which  fall  under  the  general  categories  of  approximation 
with  exponential  sums,  linear  prediction,  parameter  estimation  or  system 
identification  even  though  these  topics  are  closely  related  to  Prony's 
method.  The  papers  which  we  found  useful  in  these  fields  are  referenced 
in  the  other  volumes  which  make  up  this  report. 

Section  2.0  of  this  volume  contains  a  translation  from  the  French  of 
Prony's  original  paper.  This  has  been  included  for  two  reasons.  The 
first  being  that  while  Prony's  original  paper  is  highly  referenced,  we 
suspect  that  most  people  have  never  seen  it  let  alone  have  read  it.  The 
second  reason  for  including  it  is  that  we  find  it  extremely  interesting 
to  read. 

Section  3.0  contains  an  annotated  bibliography  of  journal  articles, 
reports  and  theses  that  deal  with  Prony's  method.  In  most  cases,  the 
annotation  was  taken  from  the  author's  abstract  of  his  psper.  Section 
4.0  contains  lists  of  papers  presented  at  conferences,  symposia  and 
meetings.  Section  5.0  contains  a  list  of  books  dealing  with  the  subject 
of  Prony's  method. 

It  is  hoped  that  the  users  of  this  list  will  cross  reference  it 
with  their  own  bibliographies  on  Prony's  method  and  will  let  ETI  know 
of  any  additions  which  they  find. 
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2.0  TRANSLATION  OF  PRONY'S  ORIGINAL  PAPER 


EXPERIMENTAL  AND  ANALYTICAL  ESSAY 
ON 

THE  EXPANSION  PROPERTIES  OF  ELASTIC  FLUIDS 
AND  ON  THE  FORCE  OF  EXPANSION  OF 
WATER  VAPOR  AND  ALCOHOL  VAPOR  AT  DIFFERENT 
TEMPERATURES 

By  R.  Prony,  1795 
Translated  by  Dr.  Ann  Sanders/ETI 

GENERAL  CONSIDERATIONS 

In  approximately  the  last  40  years ,  Physics  has  been  enriched  greatly 
by  a  large  number  of  observations  made  with  extreme  care  by  clever  and  ex¬ 
perienced  men.  This  knowledge  la  Increasing  dally,  and  Its  substance  Is 
becoming  ever  more  precious  as  tha  perfection  of  new  Instruments  provides 
more  and  more  precision  to  the  experiments;  phenomena  have  been  correlated, 
compared  and  classified;  the  language  of  an  Important  part  of  science  has 
become  analytical;  reasonable  theories  have  replaced  futile  and  often 
absurd  theories,  which  had  been  taught  until  the  middle  of  this  century. 

The  study  of  nature  in  terms  of  its  observation  and  the  knowledge  of 
Its  operation,  seems  to  me  to  consist  of  two  aspects  which  should  not  be 
confused:  the  explanation  of  effects  and  their  measurement . 

The  axplanation  of  effects  consists  of  findings,  in  a  class  of  multi¬ 
faceted  phenomena,  those  simple  or  basic  phenomena  from  which  tha  others 
are  derived,  or  of  which  they  consist  in  various  combinations.  It  also 
consists  of  showing  how,  from  the  most  varied  appearance,  one  can  untangle 
the  np«ration  and  the  being  of  the  elements  considered  to  be  the  basis 
of  the  system. 
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Thus,  by  starting  from  the  affinity  of  certain  substances  which  are 
considered  as  basic  phenomena,  it  has  been  found  that  meteorological 
phenomena,  combustion  phenomena  etc.,  result  from  these  same  affinities, 
and  show  up  under  different  disguises  which,  until  recently,  hid  their 
nature  from  physicists.  It  is  by  this  decomposition  of  complex  effects 
into  basic  effects  that  we  can  discover  certain  secrets  of  aature,  which, 
while  it  allows  us  to  lift  one  corner  of  the  veil  which  covers  it,  keeps 
another  corner  tied  by  a  knot  which  our  hand  cannot  untie. 

The  measurement  of  effects  is  the  evaluation  of  the  degree  to  which 
each  effect  is  variable  when  one  varies  either  the  causes  which  produce 
it,  or  other  effects  to  which  it  is  tied.  It  is  known,  for  instance,' 
that  the  tendency  of  fluid*  to  vaporise  is  on  the  one  hand  enhanced  by 
temperature,  and  on  the  other  hand  reduced  by  the  atmospheric  pressure, 
and  that  vaporization  occurs  only  when  the  former  force  exceeds  the 
latter  one;  in  addition  to  this  fact,  it  may  be  necessary  to  give  the 
values  of  the  pressures  which  produce  vaporization  equilibrium  at  different 
temperatures.  This  same  kind  of  reasoning  is  applicable  to  infinitely 
many  other  phenomena. 

It  can  thus  be  seen  that  the  explanation  of  effects,  whose  great 
advantage  is  to  simplify  science  and  to  correlate  its  various  aspects, 
by  analysis  and  by  decomposition  into  basic  phenomena,  must  be  complemented 
by  the  measurement  of  these  same  effects,  which  is  always  very  useful  and 
often  indispensable  when  it  is  desired  to  apply  theoretical  discoveries 
to  the  needs  of  society. 

Only  experiments  can  provide  the  initial  data  on  the  measurement  of 
physical  effects;  following  those,  however,  calculations  can  he  very  use¬ 
ful,  either  to  obtain  results  intermediate  to  those  obtained  from  the 
experiments,  or  to  correct  experimental  anomalies.  The  method  used  in 
that  case  is  called  interpolation.  Its  aim  is  to  find  a  relationship 
between  two  or  three  variables,  such  that  when  a  certain  value  is  given 
to  one  or  two  of  the  variables,  a  equally  specific  value  is  obtained  for 
the  second  or  third  variable  when  the  problem  is  considered  in  that  manner, 
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it  can  be  solved  in  an  infinite  number  of  ways,  since  there  are  an  infinite 
number  of  functions  which  can  be  obtained  by  the  same  substitutions.  It 
would  be  erroneous,  however,  to  think  that  all  those  solutions  are  equally 
applicable  to  a  given  case.  Even  though  nature  is  governed  by  seemingly 
few  and  simple  general  laws,  it  has  as  many  modifications  in  its  workings 
as  variety  in  its  forms,  and  each  measurable  phenomenon  is  always  corre- 
lated  to  a  specific  function  which  represents  it  exclusively. 

The  problem  of  Interpolation  thus  has  two  distinct  parts;  the  first 
part  consists  of  satisfying  given  numbers,  and  the  second  part  consists 
of  finding  the  function  applicable  to  the  phenomenon  amongst  all  the 
functions  which  satisfy  the  conditions. 

In  ^19  of  my  Analysis  Lessons,  I  have  given  a  solution  to  the  first 
part  of  the  problem  which  is  used  very  often,  especially  as  a  method 
of  correction.  On  the  same  subject;  Lagrange  has  published  a  beautiful 
treatise*  where  he  handled  the  subject  in  a  more  general  way  than  had 
been  dona  before.  The  students  which  know  the  theory  explained  in  my 
lessons  if  19,  20  and  21,  will  be  able  to  study  this  work  with  ease, 
and  will  be  able  to  gain  much  from  tha  time  they  will  have  spent  on  it. 

From  the  knowledge  we  have  now,  the  solution  to  the  second  part 
does  not  seem  to  be  governed  by  general  rules,  especially  whan  the  number 
of  observations  is  small,  and  whun  the  observations  do  not  span  a 
wide  range.  The  only  guidelines  available  in  this  painful  search  are 
careful  attention  to  all  the  details  and  tha  progression  of  the  experi¬ 
ments,  repetition  of  measurements,  and  analogies.  These  difficulties, 
combined  with  the  difficulty  of  precision  in  experiments,  are  the 
reason  why  exact  determinations  of  physical  laws  are  very  rare  in  Physics. 

In  1790,  1  had  the  opportunity  to  follow  very  detailed  and  careful 
experiments  on  the  expansion  force  of  water  vapor,  and  it  was  my  aim 
to  find  the  applicable  formula.  The  regularity  of  the  data  led  me  to 
believe  that  the  task  would  be  easier  than  it  actually  was.  However, 
after  some  effort,  I  found  a  certain  function  which  not  only  perfactly 
expressed  the  relationship  between  the  temperature  and  the  spring  of  the 


*  Footnote  illegible. 


2-3 


watery  vapor,  but  which  seemed  to  apply  in  general  to  phenomena  depending 
on  elastic  fluids.  I  applied  them  to  careful  experiments  which  Prieur 
made  on  the  expansion  properties  of  air  and  of  different  air-form  fluids; 
this  strengthened  my  opinion,  and  I  decided  to  publish  my  results. 

The  first  observation  which  directed  me  to  the  true  form  of  the 
function  was  the  consideration  of  certain  geometric  progressions  which 
are  present  in  certain  phenomena  related  to  elastic  fluids.  One  of  the 
most  remarkable  examples  of  this  is  the  relationship  between  the  density 
of  atmospheric  layers  and  their  respective  height.  Since  this  relation¬ 
ship  is  exponertial,  I  suspected  that  in  other  cases,  where  one  such 
quantity  would  not  be  sufficient,  two  or  more  could  be  introduced.  This 

X  X  X  X 

led  to  an  equation  of  the  form  t  ■  y,p.  +  p,p~  +  n.p,  +  ...  +  p  p 

li  44  j  j  nn 

where  t  and  n  are  the  two  variables  and  p^ ,  P^»  P ^ •••  are  the 

constants  appropriate  to  the  given  phenomenon. 

It  is  known  that  the  above  equation  results  from  an  equation  with 
linear  finite  differences  or  gives  the  general  term  in  a  recurrent  series 
of  order  n.  However,  series  of  the  type  where  a  given  term  is  obtainable 
from  a  certain  number  of  preceding  terms  seams  indeed  to  be  compatible 
with  natural  effects  where  elasticity  plays  a  large  part.  The  conserva¬ 
tion  of  living  forces  which  is  contained  in  this  property  of  bodies  always 
causes  the  actual  state  to  be  dependent  on  the  previous  states.  The  re¬ 
search  of  Lagrange,  which  I  mentioned  above,  is  also  based  on  recurrent 
series.  He  has  given  several  methods  to  find  those  which  must  interpolate 
a  given  series,  where  one  can  observe  the  elegance  and  the  depth  which 
one  comes  to  expect  from  such  a  great  analyst.  Since  the  method  I  have 
used  in  my  calculations  differs  from  his  methods,  which  I  did  not  know 
when  I  started  this  work,  I  will  now  describe  it. 
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Interpolation  Method  Applicable  to  Phenomena  Depending  on  Elastic  Fluids 


The  experiments  must  be  directed  as  much  as  possible  towards  yielding 
equidistant  results;  when  this  condition  cannot  be  met  (which  must  happen 
very  rarely) ,  but  when  the  results  are  nevertheless  numerous  and  close 
enough  together,  they  can  be  made  to  be  equidistant.  This  can  be  done 
either  by  graphic  methods,  by  drawing  the  experimental  curve,  or  by  cal¬ 
culations,  namely  by  considering  three  consecutive  results  Z  ,  Z  and 

i  xi 

Z  where  the  second  and  the  third  are  a  distance  x  and  x  respec¬ 
tively  from  the  first.  The  result  Z,  which  belongs  to  the  series  of 

equidistant  points,  and  which  is  a  distance  x  from  Z.,  can  be  calculated 

th 

from  the  following  equation,  which  is  obtained  from  those  in  my  19 
lesson  of  analysis: 


‘iii 


) 


This  formula  can  be  simplified  greatly  by  calculating  the  difference 
between  Z  and  Z1.  By  setting  Z ^  -  Z^  ■  w*  and  Z ^  ~  z±  m  oae 
obtains 


2  -  Zi  -  ~ 


x  -  x 


ii  .  i  ,, 

,x  -  x  i  x  -  x  iiv 
( - j —  ui  -  jj-  w  ) . 


In  order  to  ■**>■  Id  mistakes,  it  is  advisable  to  apply  both  the  formula 
and  the  graphic  method,  which,  when  applied  carefully  and  on  a  large 
scale,  will  give  a  precision  comparable  to  the  experiments  themselves. 


Once  this  preparation  has  been  made  (and  as  I  said,  the  cases  where 
this  is  necessary  will  be  extremely  rare) ,  a  certain  number  of  equidistant 
results  will  be  taken  which  will  encompass  either  the  entire  experimental 
range,  or  a  large  part  thereof.  Then  the  variable  Z  is  taken  which 
indicates  the  measurement  of  successive  effects  which  correspond  to  cer¬ 
tain  values  of  another  variable  x,  which  indicates  to  which  term  on  a 
given  scale  the  effects  Z  approach. 


To  satisfy  a  number  2n  of  results,  one  will  generally  have 


Vi  +  uiiBii+Vipu  + 


•  +  u(n)p(n) 


and  to  satisfy  a  number  2n  +  1  of  results, 

2  “  “lPl  *  “li»U  +  -  •  •  +  u(n)'(n>  +  %  +  1) .  (2) 

where  u^,  y^,  y ,  etc.,  and  p^,  p  ,  ®tc.,  are  constants  whose 
values  are  determined  from  the  experimental  results,  as  will  soon  be  seen. 


I  have  given  two  general  formulas,  although  usually  one  of  the  two 
would  be  sufficient,  but  1  had  a  simplification  in  mind  which  is  an 
important  inclusion  in  the  method,  and  which  will  now  be  described.  The 
determination  of  p^,  p^,  p^^,  etc.,  depends  on  the  solution  of  an 
equation,  and  by  using  the  second  formula,  one  satisfies  an  odd  number 
of  observations  by  an  equation  which  is  of  the  same  order  as  that  required 
for  satisfying  the  next  smaller  even  number  of  observations.  Thus,  one 
can  satisfy  four  or  five  observations  by  calculating  a  second  order 
equation,  six  or  seven  observations  by  calculating  a  third  order  equation, 
eight  or  nine  observations  by  calculating  a  fourth  order  equation,  and 
so  on.  It  will  almost  never  happen  that  eight  or  nine  results  must  be 
included  in  the  formula,  and  without  exceeding  the  limits  where 
methods  for  the  solution  of  numerical  equations  exist,  one  will  be  able 
to  treat  all  the  cases  which  Physics  ordinarily  presents.  In  addition 
to  that  advantage,  it  should  be  mentioned  that  there  will  be  at  most  only 
half  as  many  values  of  Z  as  there  are  observation  points,  instead  of  with 
the  formulas  which  refer  to  parabolic  curves,  and  which  always  have  as 
many  terms  as  there  are  observation  points. 


Here  then  is  the  method  by  which  the  constants  In  equations  (1) 
and  (2)  can  be  determined,  using  the  given  results. 


First  case,  where  the  number  of  observations  is  even: 


In  my  analysis  lessons  numbers  20  and  21,  I  have  shown  that  equation 
(1)  produced  the  general  term  in  a  recurring  series  of  order  n,  end  also 
that  the  terms  in  such  a  series,  taken  at  certain  equal  intervals,  always 
produced  series  of  the  same  order.  Given  that,  let  there  be  the  following 
two  series,  the  first  of  which  produces  the  observed  results  of  the 
particular  values  of  Z  furnished  by  the  experiment  and  the  second  of  which 
consists  of  the  corresponding  values  of  x,  namely  the  observed  results 


V  V  Zii*-,,Z(n);  Z(n+1) *  Z(n+2) *  ** *Z2n-l 


and  the  corresponding  values  of  x 


0;  x^;  2xi>...nxi;  (n+l)*^;  (n+2)x^. . . ,(2n-l)x^. 


The  quantities  ZQ,  Z^,  Z^,  etc.  must  form  a  recurring  series  whose  order 
must  be  found. 


Let  Aq,  A^,  a^  be  undetermined  coefficients,  satisfying  the  relation¬ 
ships 


Z0  +  A1  Z1  +  Aii  ZU  + . +  A(n)  Z(„)  ■  0 

Zl  +  A1  Zli  +  Ali  Zlil  + . +  Afn)  Z(n+1) 

Zil  +  Ai  Zlll  +  Ail  Ziv  + . +  A(n)  Z(n+2) 

Zlii  +  \  Zlv  +  Alt  Zv  + . +  Z(n)  Z(n+3) 


0 

0 

0 


An  2 .  +  A.  Z ,  .  +  A  Z ,  ,  *  + 

0  (n-1)  i  (n)  ii  (n+1) 


•+  A(n)  Z(2n-1)  ' 


(For  convenience,  A^  may  be  set  to  1  in  the  numerical  applications.) 


Aq  ^  A.  (  A,  ^ 

These  a  equations  will  produce  u  ratios  -v—  *  T~~"~  ‘  ~  *  *  *  * 

*(n)  A(n)  A(n)  A(n) 

which  make  up  the  desired  order  of  the  series,  and  one  will  obtain: 


For  n  ■  1 

*0  !  A1  ■  -h  :  zo 

Value  obtained  from  two  observations. 


For  n  “  2 


A  .  A  .  ZiZiii  ~  ZllZli 
0  '  11  20zil  -  Zi  Z1 

A  .  A  «  Z1Z11  ..l  Ylil. 

1  '  11  Vii  -  zizi 

Values  obtained  from  four  observations. 


For  n  *  3 


A_  •  A 


„  -  ZllWZ»i  +  (ZAv  -  ZilZ»l>Zlv  J-  (ZliZli  - 

111  *(Z1UZ111  *  ZliZlv)Z0  *  (Z<Z-  -  - 


t  lv  li  Hi'  i 


ii  ii 


ZiZiii,ZH 


Ai  1  Aiil 


+(ZlZiv  ~  2liZiil)Zlll  +  (ZliZll  ~  Z0Zlv)Zlv  *  (Z0Z111  ~  ZlZil l  Z' 


-(ZiZiv  '  ZiiZtii>Zi  -  (ZiiZil  -  Z0Ziv>Zii  '  <Z0Ziii  -  ZiZii>  Ziv 


.  A  +(ZHZ»  -  ZlZiii>ZiM  *  <Z0Ziii  •  ZlZll)Zlv  +  (Z1Z1  -  Z0Zli)Zv 

Au  :  A111  •  -(zliZli  -  zl2lll)zu'  -  (z0ziU  -  ziZll)zlu-  a1zl  z0zlitnlv 


Values  <•  stained  from  six  observations. 


Then,  upon  solving  the  equation. 


Ari  +  A,  a  +  A.,  a  +  A.,, a  + . .  +A,.a  -  0, 
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i  i  i  1 

the  n  roots  which  will  be  obtained  will  be  the  n  values  p.  ,  p._,  . 

1  il  ill  vn; 

from  which  the  values  of  P±,  P^,  P  ^  etc*  wi31  be  obtained.  Finally,  the 
quantities  v  ,  p  ,  p^,  etc.  will  be  given  by  the  equations: 


,  « -  p*<>  <z  -  p^>  (z  -  ^ . <j  -  p™> 

i  *,  *,  H,  \  *t  *,  *4  *4 

(°1  -  ',!>  (»i  <pi  -  . "(pl  -p(a)> 

(Z  -  p*1)  a  -  O444)  q  -  »4V>  . <2  -  pfnV 

Jii  X  X  X  X  X  X  X.  X 

(pii  "  Pi  *  *Pii  “  piii^  (pii  "  piv^ . (pii  ”  P(n)^ 

fj (2“piv)— -(z-p^> 

]iii  x  x.  x  x  x<  x  x  x 

(piii  "  pi  '  (piii  "  PiiJ  (piii  “  pivJ . (piii  “  P(n)) 


# 


I 


> 


t 


'(n) 


X  X  X 

(Z  -  P11)  (Z  -  p^)  (2  -  Pi^±) 


,(z  ~.W. 


(P 


(n) 


-  pi1>  <B(i)  - 


■’ll5  (pm  - 


pui) 


(p 


~i  _ 
(n) 


W 


by  noting  that  in  the  numerators  all  the  coefficients  of  the  powers  of  Z 

must  be  changed  to  subscripts  of  the  same  order,  i.e.,  that  Z^  must  be 

replaced  by  Z  (or  that  all  the  terms  which  do  not  contain  Z  must  be 
0  2 
multiplied  by  Zq),  that  Z  must  be  replaced  by  Z^ ,  .hat  Z  must  be  replaced 

by  Z^,  etc.,  thus  one  will  have  in  the  case  of 

n  -  l...ui  -  ZQ 

To  satisfy  two  observations. 


I 
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To  satisfy  six  observations 


This  agrees  with  the  formulas  which  I  gave  in  my  20th  lesson  of  analysis 
when  n  is  set  to  1.  In  order  to  make  t'  elimination  calculations  which 
give  u , ,  etc.,  identical  to  those  p.usented  in  the  lesson,  one  can 

Xi  Xi 

set  pi  "  »  etc.,  and  one  can  look  for  etc.,  in 

terms  of  i|>^»  etc.  One  knows  that  here  n  is  the  constant  increase 

in  x,  or  Ax.  Once  the  quantities  p^,  p^,  piii*  etc.  and  ui»  etc. 

are  determined,  their  values  can  be  substituted  the  equation  (1). 


z  •  “i  “l  *  Uli  PU  +  "m  pili  +- 


“(ti)  °(n) 


which  will  then  be  diapoead  to  satisfy  the  2n  given  observations*  and 
which  will  be  able  to  produce  some  result  intermediate  between  those 
obtained  in  fact. 

Second  case,  where  the  number  of  observations  is  odd: 

In  order  to  solve  the  second  case,  it  should  be  noted  that  aquation 

(2)  differs  only  by  the  constant  *rom  equation  (1).  Thus  the  series 

of  u  obtained  from  (2)  is  of  the  same  nature  as  that  obtained  from  (1) ,  the 

only  difference  being  that  in  (2)  each  term  is  increased  by  If 

these  terms  are  therefore  decreased  by  p,  the  remainders  will  satisfy 

(nr  l) 

the  relationships  which  make  up  a  sarles  of  order  n*.  Thus,  by  preserving 
the  notation  of  the  previous  article,  one  has 

W^n+l)5  +  Ai*Zi“u<afl)*  +  Aii^ii~W(n+l)^  +**,+  A(n)  (Z(n)“u(n-fl)J  -  0 
VW+l)5  +  Ai(Zii“1J(n-fl))  +  Aii(Ziii"M(n+l)  +,,,+  A(n)(Z(n+l)*,J(rH-l))  “  0 
VZii"M(nU))  +  Ai<Ziii“ll<n4-l))  +  Aii(Ziv“W(rrt-l))  +,,,+  A(n) (Z(iH-2)“u(nfl)>  "  0 
A0<Ziii’w(n+l))  +  Ai(Ziv“M(nfl))  +  Aii<Zv“M(rH-l))  +,,,+  A(n) (Zn+3)'M(nfl)>  "  0 


ti 


An(Z 


0N“(n-l)"ll(n+l))  +  AiCZ(n)“w(n+l))  +  Aii(Z(rH*l)“w(nfl))  +,t  ,+A(n)  (Z(2n-l)“w(n+l))“°  1 


)+  A,  (2,^-u 

,i 

'(1 

If  these  aquatlous  are  subtracted  from  one  another,  u 


A0(Z<n)"W(n+l))  +  Ai(Z(n+l))“M(n+l))  +  Aii(Zn+2))"U(n+l)>  * •+A(n) (Z(2n)"M(nfl))“° 


eliminated,  and  the  equations  will  become 


(n+1) 


will  be 


*1 


*  The  recurring  series  is,  in  fact,  of  order  n+1,  but  the  equation  of  the 
relationship  has  a  root  equal  to  unity:  l.e.,  in  the  term  y^n+1^  ^(n+l) 

one  has  p  , .  ■  1,  such  that  the  order  of  the  equation  may  be  decreased 
by  one.  n+1 
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V  20  +  Ai  1  zi  +  A11  4  zii  +""+  A(n)4Z(n)‘° 

A0  4  Z1  +  A1  1  Zli  +  AH  4  Zili  +,"-+  A(n)  4  Z(n+l)  '  0 
Ao  4  Z11  +  \  4  ZU1  +  AU  4  Zlv  +—+  A(n)  4  2(n+2)  ‘  0 
A0  4  ZU1  +  A1  4  Zlv  +  Aii  4  Zv  +—+  A(n)  4  2(nA3)  -° 


A0  4  Z(n-1)  +  Ai  4  Z(n)  +  AU  4  Z(n+1)  +,",+  A(a)  4  Z(Za-l)  "  ° 

For  cha  aaka  of  convenience  ■  on*  can  act  A  ■  1  in  tha  numarlcal 
application*. 

A~  A,  A., 

From  thoaa  equation*,  tha  valuaa  r*“  *  -f— -  ,  v  - —  ,  ate.,  can  b* 

A(n)  A(n)  A(n) 

obtainad,  namaly: 


For  n  •  1. . .  Aq  i  ■  A  i  A  Zq 


Valua  obtainad  from  thraa  obaarvationa. 


For  n  ■  2. . 


A  Z .  A  Z.  .  .  -  A  Z..  A  Z , 


a  .  A  B  ***  **  ■ 

A0  •  Ali  II0UU  ■!!,  A 


A1  1  Ali 


4  Zi  4  Z11  -  4  Z0  4  ZU1 

4  zo  4  zii  -  1  zi  1  zi 

Valuaa  obtainad  from  flva  obaarvationa. 


For  n  »  3. . . 


>(A2111AZtli-AZtlA2tv)AZli>(AZlA7.lv-A211AZui)A2ivMAzii4Zu-Az1AZl)Azv 

0'AU1  -AzillAr;ij-Ai"AZ;v)Azj-iAZ1AElv:AZltAZ1115AZ1-(AZ11AZu-AZ1AZul)A211 


2-12 


.  .  ♦  ;4Zl4Zly,~4Zil4Zlll)42lll't'(4Zll42ll~aZ0t2lv'4Zlv'l'(4Z04Zlll~4ZtaZi:>4Z' 

if  _lf  19  \  19  _/l9  19  _l9  A  9  \  A  *»  _/l9  A  V  —AW  i9  N  »9 


*u!*ui 


'(“iwiv-“iiMi1i)i2r<“ii121r4Vz1v)4Zii-(“0“iu-“i«ii,4Ziti 

^i»a-tf1alil)iitil*WiBa111-a1a11)a1,*(«1a1-tfa«11)tfa 

'(wiiiZit-4ZiiZiii)“ii-(4ZoAZii.-42i4Zii)iZiii-(“i42i-12o“u)il2iT 


Valuaa  obCAlnad  from  aavan  obaarvatiops.. 


Than,  If  k  la  a  poaitivc  numbar  amallar  chan  n,  ona  will  bw  abla  to  avaluata 
u(n4>l)  *ron  on*  c^*  n  •quationa  which  ara  comprlaad  In 


(n+1) 


JA0  2lk)  *  *1  W)  *  Au 


<A0  +  At  +  Au  + . +  A(n)>  .  A(n) 


which  muat  glva  tha  a  ana  valua  for  no  mat  tar  which  of  tha  valuaa 

of  0,  1,  2,  3...n  ona  takaa  for  k.  If  ona  than  aolvaa  tha  aquation 


‘o*V*  Au  “*  +  *iu  aS  +- 


.+  Ajn)  ■  0 


tha  n  roota  wh?.ch  will  ba  obtalnad  will  ba  tha  valuaa  of  pxi  pxi  pxl  . ..pxl 

1  *  11*  ill  (n)’ 

and  tha  valuaa  of  y^,  •tc**  oi  th*  aquation  will  ba  calculatad 

from  tha  following  foraulaa**! 

X  X  X  X 

(Z  -  p11)  (Z  -  p1j1)  (Z  -  p^)  ....  (Z  -  p(*})  (Z  -  1) 

U1  "  X.  X,  X.  X,  V  X  X.  X  X 

(Pi1  -  P^)  (Pi1  -  Piji>  ki  -  9ty)  ....(Pi1  -  P(J))  (Pi1  -  1) 


-  X  *i 

**  x$!  in  thoaa  aquatlona,  ona  dalataa  tha  factora  Z-l  and  p^  -1,  p^-i, 
p^-1,  ate.,  ona  will  hava  tha  conatanta  which  ara  appropriata  to  tha 
ganaral  tarn  of  tha  aarlaa  AZQ,  AZ^ ,  AZ^,  ate.,  i.a.,  ona  will  hava  tha 
conatanta  which  ahould  multiply  Xi  xi  ate.,  in  tha  valua  of  AZ. 

Pi  ,  PU  , 
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U(W1)  ■  2o  *  (“l  +  “ti  +  "til  + . +  %)> 

Obaarve  that,  aa  in  tha  pravioua  articla,  In  tha  coefflcianta  of  tha 

powara  of  Z  In  tha  numerator*,  auperacrlnta  of  tha  aana  number  ahould  ba 

aubatitucad,  or  that  tha  tarma  not  containing  Z  muat  ba  multipliad  by  Z.  and 

2  3  u 

thoaa  containing  Z^,  Z^,  Z  ate.,  muat  ba  multipliad  by  Z,  Z  ,  Z  ,  ate. 


Tha  laat  valua  of  ia  much  aaaiar  to  calcuiata  from  tha  pravioua  ona, 

which  ona  can  uaa  only  for  verification.  If  ona  glvaa  different  valuaa  to 


n,  ona  will  have  for 


The  above  Is  still  in  agreement  with  the  formulas  in  my  20th  lesson 
of  analysis!  when  in  each  case  m,  and  the  p  with  the  highest  superscript 
are  set  equal  to  unity. 


Once  the  numbers  p^,  p^,  p^,  •tc.,  and  y 


i*  wii*  yiii* 


etc.  are 


obtained  in  this  way,  they  can  be  inserted  into  equation  (2/ 


2  •  “t  “l1  +  “ll  #U  +  '‘ill  »iii  »lil  + . +  "(n)  p(n) 


(irt*2) 


which  will  satisfy  the  2n+l  given  observations,  and  can  be  used  to  calculate 
all  the  values  intermediate  to  those  observations. 


I  do  not  consider  the  case  where  the  equation 

A„  +  A,  a  +  A,,  a2  4- . . 4  A.  v  an  ■  0 

0  1  li  (n) 

has  jual  or  imaginary  roots;  in  my  20th  lesson  of  analysis,  I  have  given 
the  necessary  formulas  for  finding  the  solutions.  It  is  known  that  equal 
roots  Introduce  variable  and  rational  coefficients  in  the  value  of  Z,  and 
if  these  roots  are  equal  to  unity,  Z  will  contain  entirely  rational  terma: 
thus  the  interpolation  formulas  which  relate  to  rational  functions  without 
variable  dlvlsiors,  or  to  parabolic  curves,  are  only  a  very  special  case  of 
those  which  I  have  just  covered. 

I  will  now  go  on  to  applications.* 


*  The  applications,  while  interesting,  are  not  directly  relevant  to  our 
purposes  and  have  therefore  not  been  translated. 
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Auton,  J .  R. ,  "Soma  Numerical  Methods  for  Exponential  Analysis  with 
Connection  to  a  Ganaral  Identification  Schama  for  Linear  Processes" ,  M.S. 
Thesis,  University  of  Kentucky,  Lexington,  Kentucky,  1980. 

A  general  identification  scheme  for  linear  processes  is  pre¬ 
sented  and  related  to  some  well-known  exponential  analysis 
methods  such  as  Prony's  method  and  the  pencil-of-f unctions 
method.  An  elucidating  explanation  of  the  source  of  error 
(bias)  for  the  general  scheme  is  given.  Some  variation  of 
the  pancil-of-functlons  method  are  given  an  intense  exami¬ 
nation  with  some  revealing  numerical  examples.  A  new  ex¬ 
ponential  analysis  method,  called  the  adaptive  method,  is 
developed  from  the  general  scheme.  Several  numerical  ex¬ 
amples  oi  the  adaptive  method  demonstrate  its  excellent 
performance  in  the  presence  at  high  levels  of  noise.  Re¬ 
commendations  are  made  concerning  the  use  of  the  adaptive 
method  and  new  directions  for  further  research. 
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! 
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Beatty,  L.G.,  J.D.  Gaorga,  and  A.Z.  Robinson,  "Usa  of  tha  Complex  Ex** 
ponential  Expansion  as  a  Signal  Representation  for  Underwater  Acoustic 
Calibration,"  Journal  of  Acoustical  Society  of  America.  Vol.  63,  No.  6, 
pp.  1782-1794,  June  19>8. 

A  modified  Prony  method  is  developed  to  provide  a  complex 
exponential  signal  representation  for  use  with  underwater 
acoustic  calibration  waveforms.  The  approach  is  tested 
using  simulated  waveforms,  and  the  simulation  results  are 
validated  with  limited  experiments.  Results  show  the  com¬ 
plex  exponential  is  an  ideal  signal  representation  for 
underwater  acoustic  calibration  because  of  its  remarkable 
ability  to  extrapolate  calibration  waveforms  beyond  the 
actual  observation  period.  Results  are  shown  of  real  re- 
(  ciprocity  calibration  experiments  that  are  valid  down  to 

25  Hz  where  the  calibration  waveforms  are  generated  from 
a  5-ms  current  ramp  and  where  a  5-ms  observation  period 
is  used.  The  period  of  the  low-frequency  limit  is  eight 
times  the  observation  period  used  in  the  calibration. 

I  Beex,  A. A.  (Louis)  and  L.L.  Scharf,  "Recursive  Digital  Filter  Design  Via 

Covariance  Sequence  Approximation",  IEEE  Transactions  on  Acoustics. 

Speech,  and  Signal  Processing.  Vol.  ASSP-29,  No.  1,  February  l$81,p.  51- 
57. 

This  paper  proposes  a  modified  solution  to  the  general 
t  nonlinear  least-squares  approximation  problem  for  digital 

filter  design.  The  authors  choose  to  approximate  tha 
covariance  sequence  associated  with  an  ideal  filter 
because  it  arises  more  naturally  in  spectrum  matching 
problems  than  does  tha  unit  pulse  response  sequence.  The 
design  procedure  resulting  from  their  approach  requires 
I  the  rooting  of  an  nc^-order  polynomial,  tha  solution  of 

an  nth-order  system  of  linear  equations,  and  the  calcula¬ 
tion  of  a  (2n-l)  point  DFT  for  the  design  of  an  ARMA 
(2n,  2(n-l))  digital  filter. 

Bemi,  A.J.,  "Target  Identification  by  Natural  Resonance  Estimation", 

I  IEEE  Transactions  on  Aerospace  and  Electronic  Systems.  Vol.  AES-11,  No.  3, 

pp.  147-154,  March  1975. 

This  paper  presents  a  target  identification  method  based 
on  an  estimation  of  the  natural  frequencies  of  oscillation 
in  transient  radar  signatures.  The  emphasis  is  placed 
*  upon  signal  modaling  and  estimation  strategy  rather  than 

relating  resonance  locations  to  physical  structures. 

Salient  features  of  this  identification  method  ire: 

1)  target  aspect  angle  is  not  needed,  2)  multiple  targets 
of  the  same  type  can  be  illuminated  simultaneously,  and 
3)  bandpass  interrogation  pulses  can  be  used.  The  latter 
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feature  le  compatible  with  existing  radar  facilities.  The 
method  is  applied  to  some  simulated  transfer  functions, 
and  factors  affecting  estimate  accuracy  are  discussed. 

Bollig.  6.,  and  K.J.  Langenberg,  "New  Methods  of  Nondestructive  Testing  of 
Materials  and  Their  Application  Especially  in  Nuclear  Engineering". 
Presented  at  the  International  Symposium,  Taper,  Saarbrtlcken ,  Federal 
Republic  of  Germany,  17-19  September  1979  (Text  in  German). 

A  new  method  originally  developed  for  radar  purposes  is 
discussed  for  application  to  the  acouatlc  scattering 
problem.  Here  the  Helmholtz  Integral  equation  is  Laplace 
transformed  with  respect  to  time  for  the  surface  pressure 
on  a  finite  rotatlonally  symmetric  acoustically  hard 
scattering  body  in  water.  The  pertinent  algebraic  aquation 
is  obtained  by  Fourier  series  expansion  with  respect  to  the 
azimuthal  variable  and  by  applying  the  method  of  moments 
with  piecewise  parabolic  base  functions  and  collocation 
tasting.  The  singularities  of  the  resulting  Inverse  matrix 
in  the  complex  Lapicce  plane  are  computed  for  the  sphere, 
the  oblate  and  prolate  spherlod.  These  poles  are  charac¬ 
teristic  for  the  geometry  of  the  scattering  body;,  they  can 
be  used  for  classification.  Applying  Prony's  algorithm  to 
experimental  data  similar  singularities  are  extracted. 

Brlttingham,  J.N.,  E.K.  Miller,  and  J.L.  Willows,  "The  Derivation  of 
Simple  Poles  in  a  Transfer  Function  from  Real-Frequency  Information", 
Lawrence  Livermore  Laboratory,  UCRL-52050,  6  April  1976. 

A  new  technique  for  extracting  simple  poles  from  real- 
frequency  transfer-function  data  is  summarized .  The 
effects  of  varying  pole-set  parameters  are  illustrated 
with  computed  results,  and  the  technique  is  applied  to  the 
transfer  function  for  field  scattering  by  a  straight  wire. 

Brlttingham,  J.N.,  E.K.  Miller,  and  J.L.  Willows,  "The  Derivation  of 
Simple  Poles  in  a  Transfer  Function  from  Real-Frequency  Information, 

Part  2:  Results  from  Real  EM  Data",  Lawrence  Livermore  Laboratory,  UCRL- 
52118,  23  August  1976. 

A  technique  for  finding  the  simple  poles  in  a  transfer 
function  F(s)  from  real- frequency  (a»jw)  data  was 
previously  reported  in  the  reference  above  and  its 
application  demonstrated  for  a  variety  of  specific  pole 
seta.  This  procedure,  which  is  the  analog  for  frequency 
of  Prony's  technique  for  time,  is  further  tested  against 
specific  pole  sets,  and  the  results  are  described  in 
this  report. 
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A  modification  of  tha  original  procedure,  which  takes 
advantage  of  the  fact  that  ?*(s)  ■  F(s*),  is  presented. 
This  new  method  is  applied  to  electromagnetic  transfer 
functions.  The  electromagnetic  poles  thus  obtained  are 
invariant  with  respect  to  excitation  and  observable 
phenomena  (e.g.,  current,  field,  etc.)  to  which  the 
transfer  function  corresponds,  as  long  as  the  system  is 
o verde termined.  Varying  the  frequency  lntsrval  over  which 
data  is  available  also  demonstrates  the  validity  of  the 
technique. 


Buckner,  H.P.,  "Comparison  of  FFT  and  Prony  Algorithms  for  Bearing 
Estimation  of  Marrow-Band  Signals  in  a  Realistic  Ocean  Environment", 
Journal  of  Acoustical  Society  of  America.  Vol.  61,  No.  3,  pp.  756-762, 
March  1977. 


Samples  of  the  quadrature  field  {Fj}  collected  on  a 
uniform  linear  array  can  be  quickly  analysed  to  determine 
the  bearings  of  targets  by  use  of  the  FFT  or  Prony 
algorithm.  These  methods  are  compared  for  many  realiza¬ 
tions  of  a  realistic  ocean  acoustic  signal  'lus  background 
plus  noise  field.  The  familiar  Monte  Ca.'l j  mettv  1  is 
used  to  quantify  the  algorithm  performance.  It  is  shown 
that  bearing  accuracy  can  be  improved  by  forming  the 
cross-sensor  field  {Gj,}  where  Gi“<FjFj*+i>.  Theoretical 
and  computational  analysis  indicate  that  time  averaging  of 
<FjFj*+i>  not  only  reduces  the  effects  of  noise,  but  also 
reduces  Interactions  of  waves  in  the  signal  and  background 
fields  resulting  in  better  array  performance. 
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Carr,  John  W. ,  "An  Analytic  Investigation  of  Transient  Synthesis  by 
Exponentials",  Thesis,  Massachusetts  Institute  of  Technology,  1949. 

This  paper  investigates  least-mean-square  approximations 
to  transient  functions  in  the  time  domain,  offering  a 
solution  to  the  problem  of  synthesis  of  networks  to  give 
specified  transient  responses.  Various  exact  solutions 
to  this  problem  are  investigated,  none  of  which  yields  a 
solution  using  a  high-speed  electronic  digital  computer. 
Further  investigation  is  mde  of  useful  "approximations 
to  the  best  approximation."  Analysis  begins  with  the 
Prony  method  of  undetermined  coefficients  for  ordinary 
linear  differential  equations  with  constant  coefficients, 
extending  this  method  to  the  case  of  Integral  equations 
as  well.  This  enables  functions  without  continuous 
derivatives  to  be  approximated  by  a  set  of  decreasing 
exponentials,  either  real  or  complex. 

Chan,  L.C.,  and  L.  Peters,  Jr.,  "Subsurface  Electromagnetic  Mine 
Detection  and  Identification",  Ohio  State  University  ElectroScience 
Laboratory,  Technical  Report  784722-1,  November  1978. 

The  goal  of  this  research  is  to  adapt  a  commercially 
available  radar  system  developed  for  the  detection  of 
pipes  to  the  detection  of  shallow  objects  and  further  to 
establish  a  target  identification  capability.  The  target 
identification  technique  selected  for  this  purpose  (from 
among  several  being  pursued  at  the  ElectroScience 
Laboratory)  is  designated  *a  the  predictor-correlator 
scheme.  This  identification  scheme  has  already  been 
extremely  successful  in  that  a  "minelike"  target  has  been 
separated  from  a  set  of  false  targets.  The  best  designed 
and  completely  tested  system  to  date  yields  100%  identifi¬ 
cation  of  the  "minelike"  target  and  a  zero  false  alarm 
rate  of  any  of  the  false  targets.  There  remains  substan¬ 
tial  Improvements  (some  of  which  are  nearlug  completion) 
that  can  be  incorporated  in  both  ths  radar  system  and  in 
the  target  identification  techniques.  However,  the  target 
Identification  capability  discussed  in  this  report 
represents  such  a  substantial  advance  that  it  is  being 
published  in  its  present  form. 

Chan,  L.C.,  "Subsurface  Electromagnetic  Target  Characterization  and 
Identification",  Ohio  State  University  ElectroScience  Laboratory, 
Technical  Report  784722-3,  April  1979. 

A  method  for  subsurface  radar  target  characterization  and 
identification  is  described.  This  method  characterizes 
subsurface  radar  targets  by  their  complex  natural 
resonances  which  are  extracted  directly  from  their 


backscattered  time-domain  waveforms.  The  difference 
equation  coefficients  associated  with  the  complex 
resonances  are  then  used  in  the  predictor-correlator  for 
target  identification.  Both  the  characterization  and 
identification  processes  are  extensively  tested  with  real 
radar  measurements  and  found  ro  yield  practical  target 
identification  performance.  The  target  identification 
process  is  simple  and  involves  only  simple  algebraic 
operations.  Based  on  the  identification  process,  a  "first- 
generation"  microcomputer  identification  radar  system  is 
implemented  for  target  identification  in  real  time.  This 
radar  system  is  found  to  yield  practical  identification 
performance. 

Chan,  L.C.,  D.L.  Moffatt,  and  L.  Peters,  Jr.,  "A  Characterization  of 
Subsurface  Radar  Targets",  Special  Issue  on  Exploration  Geophysics, 
Proc.  of  IEEE.  July  1979. 

The  capability  of  subsurface  target  identification  at 
shallow  depths  has  been  demonstrated  using  an  electro¬ 
magnetic  video  or  baseband  pulse  radar.  Real  radar 
measurements  were  collected  for  five  targets  at  a  depth 
of  5  cm  (2  in)  in  various  ground  conditions.  These 
measurements  were  processed  for  target  characterization 
and  identification.  Identification  performance  based  on 
a  single  radar  observation  was  evaluated.  The  identifi¬ 
cation  process  requires  only  simple  algebraic  operations 
and  thus  offers  the  potential  of  real-time  on-locatlon 
identification  of  subsurface  targets.  Pole  patterns  are 
extracted  from  the  backscattered  waveforms  using  Prony's 
method. 

Cho,  K.S.  and  J.T.  Cordaro,  "Calculations  of  the  Singularity  Expansion 
Method  Parameters  from  the  Transient  Response  of  a  Thin  Wire",  AFWL 
Interaction  Note  379,  September  1979. 

The  problem  of  determining  the  singularity  expansion 
method  (SEM)  parameters  from  transient  thin-wire  data 
is  examined.  A  computer  code  is  used  to  generate  re¬ 
sponse  data.  The  SEM  parameters  are  computed  from  these 
data  using  Prony's  method.  For  noisy  dafii,  it  is  shown 
that  the  parameter  values  can  be  improved  by  signal 
averaging . 

Chuang,  C.W.  and  D.L.  Moffatt,  "Complex  Natural  Resonances  of  Radar 
Targets  Via  Prony's  Method",  Ohio  State  University  ElectroScience 
Laboratory,  Technical  Report  3424-3,  April  1975. 
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It  has  baen  shown  that  axtraction  of  complex  natural 
resonances  and  residues  from  exact  transient  responses  via 
Prony's  method  is  both  efficient  and  accurate.  However, 
in  the  art  of  radar  interrogation  the  exact  transient 
response  is  difficult  to  obtain  and  we  must  be  content  with 
an  approximate  one.  Knowing  the  complex  natural  resonances 
of  objects  is  essential  to  one  type  of  radar  target 
discrimination.  In  this  paper  we  consider  the  use  of 
Prony's  method  in  extracting  complex  natural  resonances 
from  approximate  backscattered  ramp  responses  of  radar 
targets. 

A  ramp  waveform  is  used  in  radar  target  discrimination 
because  a  good  approximation  to  the  time  domain  response 
can  be  obtained  with  only  a  few  harmonically  related 
frequency  domain  responses  via  Fourier  synthesis.  By 
applying  Prony's  method  to  this  approximate  ramp  response, 
a  few  dominant  resonances  which  are  close  to  the  origin  of 
the  complex  frequency  plane  can  be  extracted.  This  is 
quite  encouraging  because  for  the  purpose  of  radar  target 
discrimination  only  a  few  dominant  modes  are  needed. 

Objects  considered  are  prolate  spheroids,  simple  wire 
aircraft  models  and  some  realistic  aircraft  models.  The 
accuracy  of  this  method  is  found  to  be  good  when  compared 
to  the  results  obtained  by  other  techniques. 

Chuang,  C.W.  and  D.L.  Moffatt,  "Natural  Resonances  of  Radar  Targets  Via 

Prony's  Method  and  Target  Discrimination",  IEEE  Transactions  on 

Aerospace  and  Electronics  Systems,  Vol.  AES- 12,  No.  5,  pp.  583-589, 

September  1976. 

The  dominant  complex  natural  resonances  of  radar  targets 
are  obtained  via  °rrn’/’s  method  applied  to  calculated  and 
measured  backseat tertc  ramp  response  waveforms.  Subject 
targets  are  spheres,  simple  wire  models  of  straight  and 
swept  wing  aircraft,  and  realistic  models  of  modern  fighter 
aircraft.  It  is  demonstrated  ilit  when  the  backscattered 
ramp  response  waveforms  are  obtained  via  Fourier  synthesis 
of  limited  spectral  range  harmonic  scattering  data,  some 
resonance  locations  at  variance  with  those  obtained  from 
reaction  integral  equation  search  procedures  are  obtained. 

It  is  also  shown,  however,  that  the  Prony  deduced  resonances 
can  be  used  successfully  in  predictor-correlator  target  dis¬ 
crimination. 


Cordaro,  J.T.,  "A  Note  on  Representing  e  Trensient  Waveform  by  e  Finite 
Sum  of  Complex  Exponentials" ,  AFWL-TR-77-46,  Air  Force  Weapons  Laboratory, 
July  1977.  Also  Air  Force  Weapons  Laboratory  Mathematics  Notes  #46. 

The  problem  of  representing  an  observed  transient  waveform 
by  a  finite  sum  of  complex  exponentials  is  considered.  A 
least  squares  technique  that  involves  overfitting  the  model 
signal  and  then  estimating  the  correct  order  is  applied  to 
simulated  data  as  well  as  data  recorded  on  photographs.  In 
all  the  examples,  the  fitted  waveform  approximates  the 
original  waveform  wall  when  both  functions  ars  plotted  on 
the  same  graph. 

Cordaro,  J.T.,  "Comparison  of  Three  Techniques  for  Calculating  Poles  and 
Residues  from  Experimental  Data,"  AFWL  Mathematics  Notes  #61,  August  1978. 

Three  pole  and  residue  calculation  techniques  are  compered 
by  applying  the  methods  to  simulated  data.  Each  method 
works  well  when  the  data  consist  of  a  damped  sinusoid  plus 
noise.  But  when  the  data  Include  effects  from  pulser 
esynchronism  and  ground  reflection,  accurate  pole  celcu- 
lation  is  more  difficult.  For  the  methods  to  produce 
accurate  results  under  some  EMP  test  conditions,  it  is 
necessary  to  use  incident  field  data  and  information  about 
the  SEM  coupling  coefficients  and  natural  modes  of  the 
body  exposed  to  EMP. 

Cordaro,  J.T.,  "Pole  Measurements  for  ths  ATHAMAS  Pipe  Test",  AFWL 
Mathematics  Notes  #56,  August  1977. 

In  the  ATHAMAS  pipe  test,  surface  current  and  charge  density 
response  data  were  recorded  for  a  metal  cylinder  illuminated 
by  EMP.  Poles  and  residues  describing  some  of  these  data 
have  been  computed  using  an  iterative  technique.  Some  poles 
can  be  Identified  as  corresponding  to  natursl  model  of  the 
cylinder.  Others  match  the  Incident  field  poles.  Most 
poles  in  the  incident  field  are  shown  to  be  due  to  ground 
reflection. 


d«  Figueiredo,  Rui  J.P.  and  C.L.  Hu,  "Application  of  a  Frequency  Domain 
Prony  Method  to  Wide  Bandwidth  Radar  Signature  Claaelflcation,"  Department 
of  Electrical  Engineering,  Rice  University,  Houeton,  Texas,  Technical 
Report  EE  7908,  September  1979. 

Let  x  and  y  be  aignala  (i.e.,  real-valued  functione  of  time) 
of  finite  duration  end  energy.  Thie  paper  develope  a  fre¬ 
quency  domain  Prony  approach  for  interpolating,  or,  in 

general,  approximating  y(t)  by  1  a^  x(t-T^);  where 

a  *  (aj.,  ...  ,  a^)  and  £  ■  (t^ . tM)  are  real 

M-vector  parameter!.  If  x  and  y  denote  the  Fourier 
tranaforma  of  x  and  y,  than  H(u>)  ■  Y(w)/X(w)  is  a  linear 
combination  of  exponentials  with  coefficients  ai;  and 
exponents  -jiiw,  where  j  ■  i  -1  .  The  Prony  method  is 
used  to  determine  at  and  ti  from  the  samples  of  H.  These 
samples  are  obtained  irom  the  (time-domain)  samples  of  x 
and  y  by  means  of  the  DFT.  Both  the  cases  in  which  these 
samples  are  noise-free  and  contaminated  by  noise  are  con¬ 
sidered.  Special  modifications  of  the  conventional  Prony 
procedure  to  taka  into  account  the  fitting  of  complex  data 
are  presented.  These  results  are  then  extended  to  the 
case  where  y  results  from  a  linear  combination  of  several 
known  signals  xj,,  each  differently  delayed.  The  case  in 
which  y  results  from  the  convolution  of  a  known  signal  x 
with  an  unknown  impulse  response  h  is  also  considered. 

A  connection  between  Wiener's  Tauberlan  theorems  and  the 
frequency  domain  Prony  approximation  are  presented  here. 

de  Flguelredo,  Rui  J.P.  and  C.L.  Hu,  "A  Fourier-Prony  Tauberlan  Approach 
to  che  Analysis  of  a  Mixture  of  Delayed  Signals,"  Department  of  Electrical 
Engineering,  Rice  University,  Houston,  Texas,  Technical  Report  EE  7907, 
September  1979. 

A  frequency  domain  Prony  approach  is  presented  for  ex¬ 
tracting  features  of  return  signals  from  targets  illumi¬ 
nated  by  wide  bandwidth  (short  pulse)  radar.  Theoretical 
details  pertaining  to  this  approach  are  described  in  the 
above  reference.  The  features  consist  of  the  relative 
delays  and  reflection  coefficients  pertaining  to  scattering 
centers  on  the  target  representing  differently  shaped  re¬ 
gions  on  the  target  surface.  The  dimensionality  of  the 
feature  vectors  thus  constructed  is  very  low  (lese  than  tan). 
Moreover,  when  used  in  the  classification  of  targets  by  a 
nearest  neighbor  classification  strategy,  such  feature 
vectors  permit  accurate  discrimination  between  targets  that 
do  not  differ  much  in  shape;  and  also  they  are  in  a  large 
measure  insensitive  to  noise.  The  results  presented  were 
corroborated  by  computer  simulations  performed  on  the  data 
base  created  by  the  coherent  X-band  short  pulse  (0.5  nano¬ 
second)  radar  at  the  Fort  Worth  operation  Radar  Range  of 
General  Dynamics  Convair  Aerospace  Division. 
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Dudley,  D.C.,  "Fitting  Noisy  Dace  with  a  Complex  Exponential  Series.," 
Lawrence  Livermore  Laboratory,  UCRL-52242,  7  March  1977.  Alao  AFWL 
Mathematics  Notes  #31. 

The  process  of  fitting  noisy  data  with  a  complex  ex¬ 
ponential  series  is  considered.  An  analysis  of  the 
classical  Prony  method  shows  that  attempts  to  obtain 
accurate  answers  in  the  presence  of  noise  meet  with 
considerable  difficulties.  Recasting  the  problem  as 
one  of  system  identification  and  applying  the  least- 
squares  method  laads  to  biased  estimates  of  the  para¬ 
meters  in  the  characteristic  equation  for  the  system. 

The  difficulty  canters  around  the  fact  that  neither 
classical  Prony  nor  least  squares  involves  any  analysis, 
identification,  smoothing,  or  filtering  of  the  noise. 

Some  techniques  for  improving  results  are  discussed, 
and  an  iterative  generalised  laast-squares  procedure, 
which  leads  to  a  noise  filter,  is  recommended  for  further 
study. 

Dudley,  D.G. ,  "Simple  Discrete  System  Modelling  in  Hilbert  Space,"  Uni¬ 
versity  of  Arisons,  1  February  1978. 

The  design  of  a  simple  discrete  system  model  for  a 
sequence  of  data  is  given.  The  modal  is  shown  to  be 
compatible  with  approximation  in  Hilbart  space.  The 
minimum  norm  solution  is  obtained.  The  result  la  a 
structure  which  accommodates  many  of  the  results  of 
system  identification  and  linear  prediction.  The 
effects  of  error  are  briefly  msntionad. 

Dudley,  D.O.,  "Parametric  Modeling  of  Transient  Electromagnetic  Systems," 
Radio  Science.  Vol.  14,  No.  3,  pp.  387-396,  May-June  1979. 

The  modeling  of  a  transient  electromagnetic  system  as 
a  single-input,  single-output,  linear,  casual  process 
is  considered.  Two  canonical  forma  are  discussed t  one 
for  impulse  function  input  and  the  other  for  more  general 
input.  Expressions  are  derived  for  errors  in  the  em  i- 
matlon  of  the  model  parameters  caused  by  errors  in  the 
process  data.  It  is  shown  chat  Prony 's  method  is  a 
special  case  of  the  impulse  function  model  and  that 
serious  errors  occur  in  attempts  to  apply  the  method  in 
a  noisy  environment.  It  is  also  shown  that  least  squaros 
estimation  of  the  model  parameters  produces  biased  re¬ 
sults  and  that  attempts  to  overcome  this  problem  have 
led  to  nonlinear  equation  seta.  Some  recommendations 
are  made  concerning  noise  identification  and  filtering 
in  order  to  Improve  results. 
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Duhamal  P.  and  J.C.  Rault,  "Automatic  Taat  Ganaration  Tachniquaa  for 
Analog  Circuit,  tnd  Syatamas  A  Raviaw,"  1BBB  TT»nyfg^ont  9R  SiKUlil 
and  Svatama.  Vol.  CAS- 26,  No.  7,  pp.  411-440,  July  1979. 

Thi.  papar  i.  both  •  raviaw  .nd  an  aaaaaamant  of  tach- 
niqu.i  praaantly  available  for  automatic  taat  ganaration 
for  analog  ayatama  tachniquaa  for  catagoriaating  tha 
problama  o£  automatic  taating  (dafinitiona,  faulta  in 
analog  ayatama,  diffarant  type,  of  ta.ta,  m‘in 
and  diagnoaia  procaduraa) ,  charactariaation  ®nd  d*j®*jj**°n 
modaa  of  analog  ayatama,  and  tha  main  aoftwara  ingradianta 
of  automatic  taat  aquipmant  can  ba  propoaad.  Savaral  tach- 
niquaa,  raapactivaly,  procaading  from  approachaa  baaad  on 
datarminiatic  and  probabiliatic  aatimation,  taxonomical 
and  topological  analyaaii  ara  datailad.  Tachniquaa  .pacific 
to  linaar  ayatama  (aavaral  of  tham  balonging  to  tha  above 
thraa  catagoriaa)  ara  daalt  with  in  a  aaparata  jjction. 

Prony’a  mathod  ia  daacribad  aa  a  mathod  of  tranafar  function 
datarmination. 
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Fletcher,  D.L.  and  C.N.  Weygandt,  "A  Digital  Mathod  of  Transfer  Function 

Calculation,"  IEEE  Trans.  Clr.  Thaory.  Vol.  CT-18,  pp.  185-187,  January  1971. 

A  mathod  is  davalopad  that  allova  tha  diract  calculation 
of  tha  least-squares  aaeimataa  of  tha  polaa  and  raaiduaa 
of  tha  tranafar  function  of  a  linaar  ayatam  whan  tha  input 
and  output  of  tha  ayatam  la  rapraaantad  in  a  aaquantial 
numerical  format.  Tha  mathod  la  particularly  adaptabla  to 
implementation  on  a  digital  computer  and  ia  quite  efficient 
for  problaaa  of  moderate  aiaa.  Tha  mathod  preatnted  can  be 
considered  an  extension  to  Frony'a  mathod  aa  applied  to 
tranafar  function  uynthaaia. 

Frotane,  D.A.,  V.A.J.  van  Lint,  Cl.  Stettnar,  R.W.  Macgurn,  and  O.M. 

Goldatain,  "Exploding-Wira  Photon  Taating  of  Skynat  Satellite, "  IEEE 

Vol.  NS-25,  No.  6,  pp.  1349-13577^ 


Thla  paper  praaanta  a  compariaon  of  part  of  tha  data  ra- 
aulting  from  a  aariaa  of  photon,  electrical  and  calculational 
exercises  on  a  structural  modal  of  tha  Skynat  I.  Calculations 
made  for  compariaon  vieh  tha  teat  data  vara  of  two  format  coda 
simulations  and  analytic  models.  Tha  coda  calculations  vara 
tha  moat  detailed  and  included  representations  of  tha  major 
thraa-dimanaional  geometrical  features  that  vara  apparent  on 
tha  Skynat  Structure  Modal  (S8M).  For  tha  must  part,  pre¬ 
dictions  vara  made  at  high  flux,  using  a  prescribed  source 
and  a  generic  pulse.  Analytic  models  vara  used  to  understand 
tha  results  produced  by  tha  coda  calculations  and  tha  experi¬ 
ments.  Simple  analytic  circuit  models  representing  electro¬ 
magnetic  modas  of  tha  system  show  promise  of  being  reasonably 
accurate  predictive  naans  for  rapid  computation  of  tha 
electromagnetic  response. 


Gardner,  C.D.,  "Expala-Fortran  Coda  for  Exponential  Approximation  by 
Laaae  Squares,"  Lawranca  Radiation  Laboratory,  UCRL-14341,  Rav.  1, 
Univariity  of  California,  Livarmora,  March  1970. 

EXPALS  is  a  FORTRAN  program  for  obtaining  a  weighted  least- 
squares  fit  of  a  sum  of  exponentials  to  be  run  on  CDC  6600 
or  IBM  7094  computers.  This  routine  fits  a  linear  com¬ 
bination  of  real  exponential  decay  functions  by  the  Prony- 
Householder  iteration  method.  This  report  presents  the 
mathematics  and  programming  necessary  to  run  EXPALS.  Also 
Included  are  output  of  examples  to  test  the  program. 

Gavel,  D.T.,  J.V.  Candy,  and  D.L.  Lager,  "Parameter  Estimation  from  Noisy 
Transient  Electromagnetic  Measurement,"  Lawrence  Livermore  Laboratory, 
Submitted  for  Nuclear  EMP  Meeting,  Anaheim,  California,  3-7  August  1980. 

This  paper  describes  the  performance  of  three  different  signal 
processing  techniques  applied  to  parameterize  a  body  from 
noise  experimental  electromagnetic  transient  response  data. 

The  techniques  range  from  the  well-known  Prony  method  to  the 
more  sophisticated  extended  Kalman  filter  and  finally  to  the 
highly  sophisticated  maximum  likelihood  identifier.  Comparison 
is  made  of  the  performance  of  these  algorithms  and  a  discussion 
of  their  tradeoffs  is  given. 

Gerck,  E. ,  "Prony  Method  Revised,"  Applied  Optics.  Vol.  18,  No.  18, 
pp.  3075,  13  September  1979. 

The  Prony  method  described  in  its  application  to  optical 
resonators  by  Slegmau  and  Miller  has  since  than  deserved 
some  attention.  This  letter  shows  that  the  Prony  method 
for  symmetric  matrix  la  all  that  is  necessary  for  any 
scalar  resonator  calculation,  regardless  of  misalignments 
or  misflguring  of  mirrors.  This  is  accomplished  by  simple 
modifications  of  the  original  algorithm.  In  this  way  a 
simpler  and  faster  procedure  is  presented  only  by  con¬ 
sidering  the  correct  scalar  product  to  be  used.  The  pro¬ 
posed  method,  however,  could  not  possibly  apply  directly 
in  the  cases  for  which  the  polarization  of  the  beam  is 
changed  *a  in  "optical  diodes"  or  axicons. 

Goodman,  D.K. ,  "Determination  of  the  Number  of  Terms  Necessary  for  a 
Class  of  Approximation  Procedures,"  in  Data  Representation.  Edited  by 
H.S.  Andersen  and  M.R.  Osborne,  University  of  Queensland  Press,  pp. 77-93, 
1970. 


This  presentation  is  concerned  with  the  following  areas: 
(1)  the  order  determination  task  when  the  basis  functions 
are  exponentials;  (2)  a  tutorial  introduction  to  a  new 
arithmetic  designed  for  digital  computation  which  pro¬ 
duces  automatic  error  bounding;  and  (3)  a  solution,  both 
theoretical  and  practical,  to  the  problem  in  (1)  and  the 
exponential  fitting  procedure. 
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Handaraon,  T.L, ,  "Matrix  Mathoda  for  Determining  Syatam  Polaa  from 
Tranalant  Raaponaa",  Univaralty  of  Kentucky,  Elactromagnatica  Raaaarch 
Raport  30-1,  May  1980,  alao  AFWL  Mathamatica  Note  066. 

Tha  problam  traatad  la  that  of  idantifying  tha  polaa  of  a 
finita  ordar  eyatam  by  obaarving  lta  tranalant  dacay  aftar 
caaaation  of  input,  for  a  limitad  tiraa,  ualng  (poaaibly) 
multipla  obaarvation  pointa  and  axparimantal  rapatition. 

Varioua  approachaa  ara  atudiad,  having  tha  common  charac- 
tarlatlc  that  a  homoganaoua  matrix  aquation  muat  ba  aolvad. 

Savaral  tachniquaa  that  hava  baan  glvan  acant  attantion  in 
tha  lltaratura  ara  conaolidatad  into  tha  traatmant,  togathar 
with  naw  raaulta  Including  an  analytical  traatmant  of  tha 
conaaquancaa  of  aaauming  an  axcaaaivaly  high  ayatam  ordar, 
darivation  of  a  atatlatically  unbiaaad  aatimata  for  an 
lntarmadlata  paramatar  in  tha  aolutlon,  naw  thaorama  on 
arror  af facta,  a  racipa  for  affactiva  uaa  of  tha  aingular 
valua  dacompoaition,  a  naw  mathod  for  auppraaalon  of 
axtranaoua  polaa,  an  alucidating  darivation  and  axtanaion  of 
tha  mathod  of  Jain,  a  naw  form  of  tha  problam  wharain 
tha  ayatam  polaa  ara  aiganvaluaa,  and  a  atudy  of  tha 
ralatlonahlp  batwaan  various  polo  idantif ication  mathoda. 

Handaracn,  T.L,,  "Qaomatric  Mathoda  for  Datarmining  Syatam  Polaa  from 
Tranalant  Raaponaa",  to  ba  publiahad  by  IE&E  Tranaaction  on  ASSP. 

Tha  problam  traatad  la  that  of  idantifying  tha  polaa  of  a 
finita  ordar  ayatam  by  obaarvlng  lta  tranalant  daoay  aftar 
caaaation  of  input,  for  a  limitad  tima,  ualng  (poaaibly) 
multipla  obaarvation  pointa  and  axparimantal  rapatition. 

Varioua  approachaa  ara  atudiad,  having  tha  common 
charactarlatic  that  a  homogenous  matrix  aquation  must  ba 
aolvad.  Savaral  tachniquaa  that  hava  baan  given  scant 
attantion  in  tha  literature  ara  consolidated  by  a 
geometric  traatmant,  togathar  with  naw  results  including 
an  analytical  traatmant  of  tha  conaaquancaa  of  assuming  an 
excessively  high  syatam  ordar,  darivation  of  a  statistically 
unbiased  aatimata  for  an  lntarmadlata  paramatar  in  tha 
aolutlon,  naw  theorems  on  arror  affects,  a  racaipa  for 
affactiva  uaa  of  tha  elagular  value  dacompoaition,  and  a 
naw  mathod  for  auppraaalon  of  axtranaoua  polaa. 

Hiller,  J.,  "Exponential  Fitting  and  General  Tachniquaa  of  tha  Prony 
Form",  in  Data  Repraaantation.  edited  by  Andaraan  and  Oaborna,  Univaralty 
of  Quaanaland  Praaa,  pp.  62—76 ,  1970. 

Tha  intention  of  tha  flrat  part  of  th.  paper  ia  to  poaa  a 
aariaa  of  queriaa  and  to  praaant  aome  axamplaa  that  earve 
to  ahow  tha  ralavance  of  exponential  fitting.  It  ia 
damonatratad  that  tha  juatif ication  for  exponential 


fitting  is  oftsn  poor.  Nonetheless  it  is  not  always  sssy 
to  posa  tha  "corract"  fitting  problem  and  exponential 
fitting  may  ba  raasonably  invoked  by  default.  Tha  sacond 
part  of  tha  papar  discussas  tha  Prony  procadura  in  tarms 
of  commuting  oparators.  Tha  aim  of  tha  traatmant  hare  is 
to  point  out  tha  questions  that  should  ba  askad  so  that 
ganaral  oparator  claaaaa  aay  ba  handlad. 

Holt,  John  N.  and  Robyn  J.  Ant ill,  "Determining  tha  Numbar  of  Tarns  in  a 
Prony  Algorithm  Exponential  Fit,"  Mathematical  Bioseianeas.  36,  pp.  319- 
332,  1977. 

Tha  problam  of  datarmining  tha  corract  numbar  of  terms  to 
taka  in  an  axponantial  sum  approximation  to  squally  spaced 
data  is  considered.  Tha  Prony  algorithm  is  used,  together 
with  singular  value  decomposition  techniques.  A  criterion 
is  established  which  is  useful  in  estimating  tha  numbar  of 
tarms  whan  tha  data  are  subject  to  normally  distributed, 
zero  mean  noise, 

Housaholdar,  A.S.,  "On  Prony's  Method  of  Fitting  Exponential  Decay  Curves 
and  Multiple-Hit  Survival  Curves,"  Oak  Ridge  National  Laboratory,  ORNL- 
455,  Oak  Ridge,  Tennessee,  February  1950. 

This  papar  first  describes  Prony's  method.  Tharaaftar  it 
is  shown  how  one  can  procaad  to  obtain  a  valid  least 
squares  fit,  and  finally  a  criterion  for  choosing  tha 
number  of  exponentials  required  is  given.  This  is  tha 
source  of  tha  Housaholdar  orthogonalizatlon  procadura. 

Hudson,  H.G.  and  D.L.  Lager,  "Observations  on  tha  Operation  of  tha  SEMPEX 
Coda,"  Lawrence  Livermore  Laboratory,  UCID-17440,  22  September  1976. 

Tho  authors'  observations  on  tha  use  of  tha  SEMPEX  computer 
coda  are  given.  Tha  issues  of  proper  sampling  rata  and 
total  sampling  interval  are  addressed.  The  effect  that 
noise  in  tha  data  ,'s  on  tha  results  and  the  ability  of 
autocorrelation,  truncation  filtering,  and  tha  sliding- 
window  technique  to  reduce  these  affects  are  also  discussed. 


Jain,  V.K. ,  "Raprasantation  of  Sequences."  IEEE  Transactions  on  Audio  and 
Elactroacouatics.  Vol.  AU-19,  No.  3,  pp.  208-215,  September  T57TT 

In  this  papar,  cha  spaca  of  aquara  summable  raal  saquancaa 
is  considared.  Insighc  is  davalopod  into  tha  structura  of 
cha  spaca  and  a  method  for  raprasantation  of  its  elamants 
is  dascribad.  Tha  mathod  is  shovn  suitabla  for  on-lina 
ioplamantation  and  is  nonitarativa;  it  can  ba  usad  for 
wavaform  dasign,  signal  faatura  axtraction,  and  discrete- 
tima  system  identification.  Tha  applications  ara  exemplified 
by  three  simulation  studies. 

Jain,  Vijay  K.,  Tapan  K.  Sarkar,  and  Donald  D.  Weiner,  "Extension  of 
Pancil-of-Functions  Mathod  to  Ravarsa-TIme  Processing  with' First-Order 
Digital  Filters,"  Department  of  Electrical  Engineering,  Rochester 
Institute  of  Technology,  Rochester,  New  York,  Technical  Report  No.  3, 
August  1980. 

In  this  presentation,  tha  data  signal  is  processed  in 
reverse-time  by  a  cascade  of  first  order  digital  filters 
to  yield  a  family  of  information  signals.  Tha  Gram  matrix 
of  these  information  signals  is  shown  to  contain  tha 
essential  information  on  tha  poles  of  the  signal.  Tha 
entire  procedure  of  tha  application  of  pancil-of-function 
mathod  is  thus  noniterative.  Examples  presented  demon¬ 
strate  (i)  noiseworthiness  in  the  representation  problem 
when  data  ara  corrupted  by  noise,  and  (ii)  tha  effectiveness 
c>£  the  method  in  the  approximation  problem. 
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Kammler,  D.W.,  "Existence  of  Best  Approximations  by  Sums  of  Exponentials," 
Journal  of  Approximation  Theory.  Vo! .  9,  No.  1,  pp.  78-90,  September  1973. 

This  thesis  studies  the  problem  of  approximation  by  sums 
of  exponentials  and  contains  a  discussion  of  Prony's 
method . 

Kammler,  D.W.,  "Prony's  Method  for  Completely  Monotonic  Functions," 

Journal  of  Mathematical  Analysis  and  Applications.  Vol.  57,  No.  3, 
pp.  560-570,  1  March  1977. 

Prony's  method  for  constructing  an  approximating  exponential 
sum  requires  one  to  select  2n  uniformly  spaced  points,  solve 
a  difference  equation  to  find  c]_,  ...  ,  cn,  find  the  roots 
of  a  polynominal  and  take  a  log  to  find  the  exponents 
\l,  ...  ,  Xn.  This  procedure  cannot  be  used  to  produce  an 
exponential  sum  for  an  arbitrary  FeC[0,»]  but  then  one 
cannot  reasonably  expect  to  model  an  arbitrary  continuous 
function  with  such  an  exponential  sum.  This  paper  proceeds 
under  the  additional  hypothesis  that  F  is  completely  nono¬ 
tonic  on  [0,-].  It  is  shown  that  Prony's  method  always 
yields  an  exponential  sum  which  for  sufficiently  large  n 
can  be  made  arbitrarily  cose  to  F  (in  the  uniform  norm)  if 
and  only  if  F  is  completely  monotonic  on  [0,“1. 

Kelso,  J.M.,  "Measuring  the  Vertical  Angles  of  Arrival  of  HF  Skywave 
Signals  with  Multiple  Modes,"  Radio  Science.  Vol.  7,  No.  2,  pp.  245-250, 
February  1972.  '  ” 

The  measurement  of  the  vertical  angle  of  arrival  of  HF 
skywave,  multimode  signals  generally  requires  either 
highly  complicated  mathematical  manipulations  or  observing 
techniques  that  necessitate  a  waiting  period  until  the 
relative  phases  among  the  modes  either  reach  a  desired  con¬ 
dition  or  pass  through  a  desired  range  of  variations.  This 
paper  utilizes  an  extension  of  a  very  old  method  (the  Prony 
method,  c.  1795)  for  approximating  sums  of  exponential  func¬ 
tions.  The  technique  uses  relatively  simple  calculations  that 
are  based  on  <axi  instantaneous  (or  nearly  so)  set  of  samples 
from  a  uniformly  spaced,  vertically  disposed  array  of  antenna 
elements.  Unlike  many  other  methods,  this  technique  requires 
neither  time  variations  of  phase  nor  a  perfectly  reflecting 
ground  plane,  and  is  valid  for  arbitrary  polarization  of  the 
(similar)  antenna  elements  and  for  any  general  wave  polariza¬ 
tion.  Experimental  results  from  some  field  tests  of  the  pro¬ 
cess  for  vertical  angle  estimation  are  presented. 

Khatwanl,  K.J.,  and  J.S.  Bajwa,  "Identification  of  Linear  Time- Invariant 
Systems  Using  Exponential  Signals,"  IEEE  Transactions  on  Automatic  Con¬ 
trol.  pp.  146-148,  February  1975. 
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A  new  method  of  system  identification  for  linear  time 
invariant  systems  based  on  Prony's  method  of  exponential 
interpolation  is  suggested  in  this  correspondence.  This 
method  of  Identification  uses  an  exponential  signal  for 
system  excitation  and  has  the  advantages  of  simplicity, 
independent  determination  of  numerator  and  denominator 
of  the  system  transfer  without  any  difficulty.  The 
simplicity  of  the  method  has  resulted  because  of  the 
simple  but  particularly  suitable  formulation  in  the  state 
variable  form  used. 

Kremer,  H.  and  A.  Lacroix,  "Praktische  Berer.hnung  von  Exponential- 
approximationen  mlt  dam  Prony-Verf ahren" ,  Angewandte  Informatik.  Vol.  3, 
pp.  121-126,  1972. 

Prony's  method  is  described  for  the  approximation  of  discrete 
sets  of  points  by  exponential  functions  with  complex  argument.1 
The  determination  of  the  degree  of  the  approximating  function 
is  discussed  in  detail.  It  is  shown  by  examples  how  input 
».rors  Influence  the  results  of  the  method.  The  algorithm 
has  been  formulated  in  ALGOL  and  tested  on  an  IBM  7040  com¬ 
puter;  the  ALGOL  program  is  given. 

Kremer,  H. ,  "Design  of  Linear  Networks  with  Prescribed  Impulse  Response", 
Extended  paper  from  IEEE  Conference  on  Computer-Aided  Design  of  Electronic 
and  Microwave  Circuits  and  Systems,  Hull,  United  Kingdom,  12-14  July  1977. 

A  method  for  the  design  of  linear  networks  with  prescribed 
impulse  response  is  presented,  that  uses  Prony's  algorithm 
for  the  approximation  And  gives  the  network  as  a  parallel 
structure  of  at  most  second  order  subnetworks.  Following 
a  brief  description  of  the  algorithm  some  numerical  end 
computational  aspects  are  discussed.  The  method  is 
illustrated  by  two  examples,  the  design  of  a  Hilbert  filter 
and  an  Ideal  low-pass  filter. 

Kulp,  R.C.,  "An  Optimum  Sampling  Procedure  for  Use  with  the  Prony  Method", 
IEEE  Trans,  on  Electomagngtic  Compatibility.  Vol.  EMC-23,  No.  2,  May  1981, 
pp.  67-71. 


The  ill-posed  nature  of  identifying  the  poles  of  a 
transient  response  is  discussed  along  with  a  sampling 
scheme  for  minimizing  the  effect  of  this  111-posedness. 

Two  simulation  examples  are  given.  The  results  of  the 
second  example  are  compared  to  previously  published 
results. 

Kunz,  K.S.,  et  al.,  "Surface  Current  Injection  Techniques:  A  Theoretical 
Investigation",  IEEE  Transactions  an  Nuclear  Science,  Vol.  NS-25,  No.  6, 
pp.  1422-1427,  December  1978. 


This  theoretical  study  was  undertaken  to  find  surface 
current  injection  technique  configurations  capable  of  quality 
simulation.  A  configuration  capable  of  quality  simulation 
was  found.  It  compared  favorably  with  conventional  EMP 
simulation  techniques  when  both  were  measured  against 
theoretical  threat  aircraft  responses.  Some  preliminary 
experimental  results  further  confirm  the  capabilities  of 
this  configuration  while  validating  the  theoretical  predic¬ 
tions.  Prony's  method  was  used  to  find  the  poles  and  residues 
of  the  data. 

Kunz,  K.S.  and  J.F.  Prewitt,  "Practical  Limitations  to  a  Natural  Mode 

Characterization  of  Electromagnetic  Transient  Response  Measurements," 

IEEE  Transactions  on  Antennas  and  Propagation.  Vol.  AP-28,  No.  4,  pp. 

577,  July  1980. 

Practical  limitations  arising  from  limited  signal-to-noise 
ratios  in  actual  transient  electromagnetic  measurements 
limits  the  number  of  natural  modes  which  can  be  extracted 
from  such  measurements.  This  number  is  shown  here  to  be 
three  to  six  for  a  complex  scatterer,  in  this  case  an 
aircraft,  five  for  a  "fat"  cylinder  and  thin  wire.  This 
is  fcr  double  exponential  waveform  excitation  of  the  air¬ 
craft  and  step  function  excitation  of  the  fat  cylirder  and 
thin  wire.  A  few  more  modes  can  be  resolved  if  the  excita¬ 
tion  is  altered  so  as  to  provide  more  high  frequency  energy, 
as  can  be  accomplished  by  changing  the  step  excitation  to  a 
delta  function  excitation,  or  by  signal  conditioning,  for 
example  making  a  derivative  measurement  that  emphasizes  the 
high  frequency  response.  Even  so,  far  fewer  modes  can  be 
resolved  experimentally  than  can  be  predicted.  Experimental 
electromagnetic  transient  response  measurements  should 
therefore  be  made  with  these  practical  limitations  in  mind. 


575- 
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Lager,  D.L. ,  H.G.  Hudson,  and  A.J.  Poggio,  "User's  Manual  for  SEMPEX:  A 
Computer  Code  for  Extracting  Complex  Exponentials  from  a  Time  Waveform" 
AFWL  Mathematics  Note  // 45,  March  1977. 

This  report  Is  the  user's  manual  for  the  SEMPEX  computer 
code.  Singularity  Expansion  Method  Pole  Extraction  Is  a 
technique  In  electromagnetics  In  which  the  free  response 
of  a  structure  is  expressed  as  a  weighted  sum  of  complex 
exponentials  where  the  damping  factors  of  the  exponentials 
are  derived  from  the  poles  (l.e. ,  singularities)  of  the 
transfer  function  of  the  object. 

Lager,  D.L.,  "Using  Parer  (Parameter  Estimation  Technique)  for  the  De¬ 
termination  of  Mode  Shapes,"  Lawrence  Livermore  Laboratory,  UCID-17898, 

8  September  1978. 

This  report  applies  the  Parameter  Estimation  Technique 
(PARET)  to  the  measured  responses  of  a  T-shaped  aluminum 
plate.  The  frequency  and  damping  for  the  natural  resonances 
and  display  animated  mode  shapes  for  the  plate  are  obtained. 

The  measurements  consist  of  the  response  of  an  accelerometer 
mounted  on  a  corner  of  the  plate  to  hammer  blows  at  38 
separate  locations  on  the  plate.  PARET  is  used  to  compute 
the  natural  resonances  by  determining  the  parameters  of  an 
exponential  model  that  fits  the  responses.  Tabulation  of 
the  frequencies  and  damping  for  seven  modes  occurring  over 
a  frequency  range  500-1650  Hz  and  also  plot  the  shape  for 
each  mode  are  made. 

Lager,  D.L.,  "A  Statistical  Evaluation  of  the  Performance  of  the  Time- 
Domain  Paret  Algorithm,"  Lawrence  Livermore  Laboratory,  UCID-18120, 

1  June  1979. 

A  methodology  was  developed  for  evaluating  the  performance 
of  the  time-domain  PARET  (TDP)  algorithm  when  processing 
waveforms  contaminated  with  additive  noise.  The  methodology 
is  general;  it  may  be  applied  to  evaluating  the  performance 
of  other  algorithms  and  provides  a  good  basis  for  comparison. 

The  performance  is  evaluated  by  computing  the  bias,  variability, 
and  root  mean  square  error  in  the  parameters  estimated  by  the 
algorithm  when  processing  ensembles  of  waveforms.  The  wave¬ 
forms  used  in  this  study  were  the  displacements  as  a  function 
of  time  from  a  computer  simulation  (by  the  SAP  IV  Program) 
of  an  eight-floor  structure. 


Lai,  Hang  H.  and  John  D.  Dyson,  "The  Determination  of  the  Direction  of 
Arrival  cf  an  Interference  Field",  Radio  Science,  Vol.  16,  No.  3, 
pp.  365-376,  May-June  1981. 
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The  problem  of  resolving  Che  complex  incerference  field  ac 
a  high  receiving  site,  that  Is  caused  by  multimode  propaga¬ 
tion  at  high  frequencies,  into  the  individual  plane  wave 
components  and  determining  the  direction  of  arrival  of 
these  components  is  considered.  The  method  developed  is 
based  on  Prony's  algorithm.  To  illustrate  the  effect  of 
noise  on  the  solution,  a  first-order  linearized  model  of  a 
statistical  study  of  a  two  plane  wave  problem  for  this 
conventional  Prony's  method  is  presented.  A  modification 
of  Prony's  technique,  an  eigenvalue  approach,  allows  one 
to  analytically  determine  the  number  of  waves  in  the 
problem,  and  this  approach  is  less  susceptible  to  the 
effect  of  noise.  Data  obtained  from  a  controlled  experi¬ 
ment  when  a  signal  from  a  distant  source  was  repeated  and 
retransmitted  at  a  different  location  verify  that  the 
system  as  well  as  the  algorithm  is  capable  of  resolving 
waves  coming  from  different  directions.  Experimental 
evidence  related  to  multimode  as  well  as  single-mode 
propagation  coming  from  the  same  source  has  been  obtained. 

Leveckis,  A.S.,  "Application  of  Singularity  Expansion  to  Electromagnetic 

Pulse  Response  Waveforms,"  Master's  Thesis,  Massachusetts  Institute  of 

Technology,  June  1977. 

Here  is  presented  the  results  of  an  attempt  to  use  Prony's 
method  on  EMP  response  waveforms  on  ships.  The  technique 
to  acquire  the  experimental  data  is  detailed;  the  numerical 
problem  is  posed,  and  the  goals  are  quickly  simplified  to 
attempting  to  find  the  best  curve-fit  to  the  waveform,  using 
complex  exponentials,  and  then  hoping  to  eliminate  the  curve- 
fit  poles.  A  method  is  presented  which  aids  in  the  'physical 
pole'  selection  process,  and  explained  also  is  the  signi¬ 
ficance  of  the  right  half  plane  poles  in  the  results. 

Lytle,  R.J.  and  D.L.  Lager,  "Using  the  Natural-Frequency  Concept  in  Re¬ 
mote  Probing  of  the  Earth,"  Radio  Science.  Vol.  11,  No.  3,  pp.  199-209, 

March  1976. 

A  technique  is  given  for  determining  the  thicknesses  and 
the  electrical  constitutive  parameters  of  a  planar-layered 
medium  such  as  .i  coal  seam  in  a  mine  environment.  Time- 
domain  experimental  data  are  analyzed  with  Prony's  method 
to  determine  the  natural  frequencies  of  the  layered  medium. 
Explicit  relations  are  given  (for  dielectric  layers)  for 
determining  the  thicknesses  and  dielectric  constants  from 
the  experimentally  determined  natural  frequency  results. 

Explicit  expressions  are  also  given  (for  conductive  layers) 
for  calculating  the  electrical  thickness  from  natural  fre¬ 
quency  results.  The  method  is  illustrated  with  sample 
numerical  results. 
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Malthoul,  J. ,  '’Linear  Prediction:  A  Tutorial  Review",  Proceedings  of  the 
IEEE.  Vol.  63,  No.  4,  pp.  561-580,  April  1975. 

This  paper  gives  an  exposition  of  linear  prediction  in  the 
analysis  of  discrete  signals.  The  signal  is  modeled  as  a 
linear  combination  of  its  past  values  and  present  and  past 
values  of  a  hypothetical  input  to  a  syster  whose  output  is 
the  given  signal.  In  the  frequency  domain,  this  is  equiva¬ 
lent  to  modeling  the  signal  spectrum  by  a  pole-zero  spectrum. 

The  major  part  of  the  paper  is  devoted  to  all-pole  models. 

The  model  parameters  are  obtained  by  a  least  squares  analysis 
in  the  time  domain.  Two  methods  result,  depending  on  whether 
the  signal  is  assumed  to  be  stationary  or  nonstationary.  The 
same  results  are  then  derived  in  the  frequency  domain.  The 
resulting  spectral  matching  formulation  allows  for  the  modeling 
of  selected  portions  of  a  spectrum,  for  arbitrary  spectral 
shaping  in  the  frequency  domain,  and  for  the  modeling  of 
continuous  as  well  as  discrete  spectra.  This  also  leads  to 
a  discussion  of  the  advantages  and  disadvantages  of  the 
least  squares  error  criterion.  A  spectral  interpretation  is 
given  to  the  normalized  minimum  prediction  error.  Applica¬ 
tions  of  the  normalized  error  are  given,  including  the  de¬ 
termination  of  an  "optimal"  number  of  poles.  The  use  of 
linear  prediction  in  data  compression  is  reviewed.  For 
purposes  of  transmission,  particular  attention  is  given  to 
the  quantization  and  encoding  of  the  reflection  (or  partial 
correlation)  coefficients.  Finally,  a  brief  Introduction 
to  pole-zero  modeling  is  given. 

l-iaria,  G.A.  and  M.M.  Fahmy,  "l“  Approximation  by  Exponentials",  IEEE 
Transactions  on  Circuit  Theory,  pp.  71-74,  January  1973. 

An  interative  method  is  proposed  for  obtaining  a  solution 
to  the  problem  of  lm  approximation  by  exponentials.  The 
method  always  converges  and  has  the  advantage  that  the 
samples  need  not  be  equidistant.  Using  the  proposed  method 
and  allowing  p  to  increase  yields,  in  a  much  easier  way, 
results  that  are  comparable  to  those  obtained  by  Fischi  for 
the  Chebyshev  approximation. 

Marzollo,  A.,  "On  the  Mean  Square  Approximation  of  a  Function  with  a 
Linear  Combination  of  Exponentials",  Int.  Journal  Control.  Vol.  9,  No.  1, 
pp.  17-26,  1969. 

This  paper  deals  with  the  problem  of  finding  the  2N 
coefficients  of  a  linear  combination  of  N  exponentials 
which  minimizes  its  mean  square  error  from  a  given 
function  f(t).  The  real  difficulty  lies  in  the  determina¬ 
tion  of  the  optimal  exponents. 
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A  real  function  of  N  real  variables  is  constructed,  such  that 
the  points  where  it  is  minimal  give  the  optimal  exponents. 

Such  a  function  may  be  considered  as  the  sample  for  t  T  of 
the  response  to  f(T  t)  of  a  simple  system.  This  implementa¬ 
tion  allows  a  significant  simplification  of  the  given  problem 
from  a  computational  point  of  view,  for  Instance  by  hybrid 
methods. 

Millar,  E.K. ,  at  al.,  "Radar  Target  Classification  Using  Tamporal-Mode 
Analysis",  Lawrence  Livermore  Laboratory,  UCRL-5182S,  27  May  1973. 

Exploratory  test  results  are  presented  for  a  new  radar 
target-classification  technique  based  on  processing  the 
transient  scattered  waveform  from  a  target  using  temporal¬ 
mode  analysis.  The  new  method  offers  the  potential 
advantages  of  target  characterization  independent  of 
aspect  angle  and  polarization,  a  means  for  quickly 
comparing  targets  against  a  stored  library,  and 
substantially  simplifying  data  storage  and  processing 
requirements.  Although  several  theoretical  questions 
remain  to  be  answered,  initial  test  results  look  quite 
promising  for  this  new  approach. 

Miller,  E.K. ,  J.N.  Brittingham,  and  J.L.  Willows,  "The  Derivation  of 
Simple  Poles  in  a  Transfer  Function  from  Real  Frequency  Information 
Part  3:  Object  Classification  and  Identification",  Lawrence  Livermore 
Laboratory,  UCRL-52211,  10  January  1977. 

This  report  is  the  third  in  a  series  of  three  that 
evaluate  a  technique  (frequency-domain  Prony)  for 
obtaining  the  poles  of  a  transfer  function.  The  main 
objective  waa  to  assess  the  feasibility  of  classifying  or 
identifying  ship-like  targets  by  using  pole  sets  derived 
from  frequency-domain  data.  A  predictor-correlator 
procedure  for  using  spectral  data  and  library  pole  sets 
for  this  purpose  was  developed.  Also  studied  was  an 
iterative  method  for  reducing  the  effects  of  noise  and  a 
technique  based  upon  magnitude-only  spectral  data. 

Miller,  E.K. ,  "Data  Characterization  and  Compression",  Lawrence  Livermore 
Laboratory,  UCID-17511,  6  July  1977. 

The  information-transformation  process  based  upon  Prony's 
method  has  been  found  to  be  increasingly  useful  for  appli¬ 
cation  to  electromagnetic-type  problems  in  particular  and 
a  growing  variety  of  physical  problems  in  general. 

Miller,  E.K.  and  D.L.  Lager,  "Information  Extraction  Using  Prony's  Method" 
Lawrence  Livermore  Laboratory,  UCRL-52329,  23  August  1977. 


This  rtport  considers  chs  questions  of  information  contsnt 
and  data  transformation.  To  dafina  information  content, 
tha  concapts  of  data  rank  (numbar  of  polas  in  tha  axponantial 
data)  and  data  pracision  ara  uaad.  To  study  data  trans¬ 
formation,  a  wavaform  is  dafinad  in  barms  of  spacifiad 
polas  and  residues,  is  samplad  to  obtain  tha  associatad 
data,  and  is  oparatsd  on  using  Prony's  mathod  to  transform 
tha  data  and  compara  it  with  tha  original  information 
contant.  Othar  numarlcal  procassaa,  such  as  matrix 
invarslon,  can  ba  vlawad  in  tha  sama  way.  Tha  problams  of 
ill-conditioning  ara  ralatad  mors  to  tha  data  than  to  tha 
procaas.  Finally,  applications  of  Prony's  mathod  to  othar 
physical  problams  with  axponantial  solutions  arc  considarad. 

Millar,  E.K.  and  O.L.  Lagar,  "Prony's  Mathod  for  tha  Angla  Domain", 
Lawranca  Livarmora  Laboratory,  UCID-17502,  Rav.  1,  6  Octobar  1977. 

Prony's  mathod,  a  tachnlqua  for  obbaining  tha  paramatars 
of  a  sarias  of  complax  axponantlals ,  is  uaad  hara  to 
analyza  tha  radiation  pattam  of  a  linaar  array  of  iso¬ 
tropic  radiators.  Tha  locations  and  amplitudes  of  tha 
radiators  ara  darivabla  from  tha  pattern.  Possible 
applications  include  imaging  and  array  synthesis. 

Millar,  E.K.  and  D.L.  Lagar,  "Radiation  Flald  Analysis  and  Synthesis  Using 
Prony's  Mathod",  Lawranca  Livarmora  Laboratory,  UCRL-B0532,  16  January 
1978. 


Prony's  method,  a  tachnlqua  for  extracting  tha  paramatars 
of  an  axponantial  sarias  from  its  dlscrata  samples,  is 
applied  hara  to  far-flald  pattern  data.  Tha  mathod  offers 
possibilities  in  imaging  source  distributions  and  in  array 
analysis  and  synthesis. 

Millar,  E.K. ,  "Prony's  Mathod  Revisited",  Lawranca  Livarmora  Laboratory, 

UCRL-52590,  18  Octobar  19""' . 

This  report  contlnuas  evaluation  of  tha  moving-win  >• 
approach  for  reducing  tha  deleterious  affects  t  f  ><>..»a  to 
which  Prony  processing  is  particularly  sensitive,  a 
sarias  of  systematic  computer  experiments  have  bean 
conducted  on  analytically  spacifiad  transient  waveforms 
and,  in  particular,  have  examined  transmission  of  errors 
though  the  process  by  determining  tha  accuracy  of  tha 
input  data,  the  characteristic-aquation  coefficients  and 
roots,  and  tha  polas.  Statistical  distributions  of  thasa 
quantities  for  combinations  of  input  data  paramatars  have 
bean  derived. 


E.K.  and  D.L.  Leger,  "Inversion  of  One-Dimensional  Scattaring 
Data  Using  Prony'a  Mathod",  Lawrence  Livermore  Laboratory,  UCRL-52667, 

12  February  1979. 

A  one- dimensional  configuration  la  tha  a imp last  gaooatry 
to  lnvart,  yat  it  haa  practical  application  to  such 
problems  aa  acattaring  from  Inhomogeneous  half-spaces  and 
propagation  on  non-uniform  transmission  lines.  Whether 
the  medium  parameters  vary  continuously  or  discretely  with 
position,  the  problem's  numerical  description  can  usually 
be  developed  in  finite- time  epproximatlon.  As  such,  the 
scattered  and  transmitted  fields  can  be  represented  aa 
exponential  series,  whose  exponents  are  related  to  the 
electrical  thicknesses  of  the  layers  which  make  up  the 
model.  If  the  exponents  or  poles  are  derivable  from  field 
data,  then  the  inverse  problem  is  formally  solvable.  Thin 
report  considers  application  of  Prony’a  method,  a  procedure  , 
for  obtaining  the  poles  of  exponential  signals,  to  such 
one-dimensional  problems.  Analysis  of  both  time-domain  and 
frequency-domain  data  is  studied.  The  effects  of  the 
medium  characteristics,  number  of  layers,  and  other  factors 
are  examined.  It  is  concluded  that  Prony'a  method  has 
merit  for  certain  classes  of  one-dimensional  inverse 
problems. 

Miller,  E.K. ,  "A  Study  of  Target  Identification  Using  Electromagnetic 
Poles",  Lawrence  Livermore  Laboratory,  UCRL-32685,  13  March  1979. 

Three  techniques  to  identify  radar  targets  by  using  their 
electromagnetic  resonances  or  poles  are  studied.  A 
singular  advantage  of  such  e  pole-based  approech  is  that 
it  can  operate  Independently  of  the  target  excitation.  A 
library  of  pole  sets  is  used  both  to  generate  the  target 
waveforms  (to  which  noise  is  added)  and  to  correlate  with 
parameters  extracted  from  those  waveforms.  The  techniques 
studied  are  linear  prediction,  residue  calculetlon,  and 
pole  calculation.  Direct  waveform  correlation  was  also 
evaluated  for  comparison  with  the  pole-based  techniques. 

Using  10- run- average  correlation  values,  error-free 
identification  was  accomplished  by  pole  calculation  at  a 
peak-slgnal-to-peak-noise  ratio  as  low  as  10,  and  by  wave¬ 
form  correlation  at  a  noise  ratio  of  unity. 

Miller,  E.K. ,  D.L.  Lager  and  J.T.  Okada,  "Imaging  of  Linear  Source 
Distributions",  Lawrence  Livermore  Laboratory,  UCRL-52822,  November  1979. 

An  approech  based  on  Prony's  method  is  used  to  find  the 
spatial  poles  of  straight  wires  of  varying  length  and 
radius  from  the  far-fleld  patterns.  For  the  antenna  case, 
poles  having  tha  largest  residues  are  found  near  the  feed 
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region  and  antenna  ends;  amallar  ones  are  found  outside 
these  regions.  The  pole*  are  attributed  to  radiation  from 
the  exciting  source  and  wire  ends,  due  to  charge  accelera¬ 
tion,  and  from  a  slightly  damped  traveling  wave.  In  the 
scattering  case,  the  poles  with  the  largest  residues  are 
located  at  the  wire's  ends  for  near-axial  Incidence  and 
are  again  dua  to  end  radiation.  By  examining  the  residue 
magnitudes  as  a  function  of  wire  radius,  a  correlation  is 
found  with  radiation  from  the  feed  region  of  an  infinite 
antenna  and  from  the  ends  of  a  semi-infinite  antenna. 

Millar,  G.,  "Representation  and  Analysis  of  Signals;  Part  XXVI,  Least- 
Squares  Approximation  of  Functions  by  Exponentials",  Department  of 
Electrical  Engineering,  The  John  Hopkins  University,  Baltimore,  Maryland, 
June  1969. 

The  approximation  of  an  analytic  time  function  in  the  least- 
squares  sense  by  sums  of  exponentials  is  considered  from 
several  different  points  of  view.  McDonough's  method  for 
finding  the  poles  is  derived  by  three  different  Approaches. 

In  the  process,  a  new  method  is  developed  which  offers  the 
advantages  of  the  earlier  results  achieved  by  McDonough  and 
by  McBride,  Schaefgen,  and  Steiglltx.  This  new  method 
reveals  the  link  between  these  earlier  methoda  and  provides 
a  standard  for  comparing  chase  two  linear  iterative  schemes 
using  several  numerical  examples. 

Mlttra,  R.  and  C.O.  Stearns,  "Field-Strength  Measurements  in  a  Multipath 
Field  Using  Linear  and  Circular  Proving",  Radio  Science.  Vol.  2,  Ho.  1, 
pp.  101-110,  January  1967. 

Techniques  for  resolving  the  components  of  a  multipath  field 
are  investigated  in  this  paper.  Both  linear  and  circular 
probing  schemas  are  discussed.  A  method  for  enalyzing  the 
linearly  probed  data  is  given,  which  is  considerably  simpler 
than  the  procedure  followed  by  Wattarson  in  an  «*»..llar  paper. 

The  method  is  alio  extended  to  the  circular  case.  It  is 
found  that  the  resolution  in  the  circular  case  is  poorer  than 
for  the  linear  scheme. 

Moffett,  D.L.  and  K.A.  Shubert,  "Natural  Resonances  Via  Rational 
Approximations",  IEEE  Transactions  on  Antennas  and  Propagation.  Vol.  AP-25, 
No.  5,  pp.  657-660,  September  1977. 

An  approech  for  locating  the  dominant  complex  natural 
reeonancaa  (poles)  of  scatt-rers  or  radiators  is  described. 

Admittance  values,  calibrated  uelng  Richmond's  thin  wire 
comouter  programs,  at  a  regular  grid  of  complex  and  real 
frequencies  ere  fit  by  a  rational  function.  Zeros  of  the 
resultant  numerator  and  denominator  polynomials  than  yield, 


respectively,  the  zeros  and  polaa  of  the  structure.  The 
method  is  illustrated  for  resonances  of  straight  end  bent 
wire  structures. 

Moffett,  D.L.,  L.C.  Chan  and  G.A.  Huwishar,  "Characterization  of 
Subsurface  Electromagnetic  Soundings",  Ohio  State  ElectroScience 
Laboratory,  Annual  Technical  Report  4490-1,  September  1977. 

In  this  report,  details  on  a  first  year's  efforts  to 
characterize  subsurface  radar  target  signatures  arw 
summarized.  The  report  covers  three  main  topics;  an 
initial  measurement  program  to  obtain  transient  radar 
signatures  of  burled  objects,  a  description  of  a  compre¬ 
hensive  computer  program  to  extract,  by  various  means,  the 
complex  natural  resonances  of  burled  objects,  end  the 
complications  Involved  in  extending  a  synthetic  radar 
Imagery  technique  to  buried  targets.  In  a  conclusions 
section,  the  state  of  the  art  is  summarized  and  research 
for  the  next  year's  effort  is  discussed. 

Moffett,  D.L.  and  L.C.  Chan,  "Characterization  of  Subsurface  Soundings", 
Ohio  State  ElectroScience  Laboratory,  Final  Report  4490-2,  January  1979. 

Research  on  a  National  Science  Foundation  Grant  for  the 
period  1  June  1976  to  31  December  1978  is  summarized. 

Unique  new  excitation  invariant  parameters  (complex 
natural  resonances)  of  subsurface  anomalies  are  described. 

Methods  for  extracting  and  exploiting  these  excitation 
invariant  parameters  in  the  presence  of  noise  and  cluttnr 
are  given.  State-of-the-art  subsurface  radar  systems  are 
illustrated.  Using  such  a  system,  target  characterization 
procedures  are  successfully  illustrated  for  a  variety  of 
subsurface  targets  using  real  radar  data. 

Moffett,  D.L.  ft  al,  "Trans Lent  Response  Characteristics  in  Identification 
end  Imaging",  IEEE  Transactions  on  Antennas  and  Propagation.  Vol.  AP-29 , 
No.  2,  pp.  192-205,  March  1981. 

Three  dlfferuut  solutions  to  problems  in  identification  and 
inverse  scattering  using  restricted  far-zons  scattering 
date  ere  reviewed  and  illustrated.  The  data  ere  restricted 
in  the  sense  that  the  frequencies,  aspects,  end  polarizations 
are  limited.  All  throe  solutions  are  based  in  part  on  a 
fundamental  time  domain  viewpoint  whereby  the  far-rone 
scattering  charactevlr tics  of  any  finita  object  are  uniquely 
summarized  by  its  impulse  response  waveforms.  Insight 
provided  by  these  waveforms  is  sxplolted  via  gsometrical 
characteristics  extracted  from  steady-state,  forced,  *nd 
free  responses  of  the  object.  For  the  present  purposes  the 
fundamental  importance  of  signaling  wave forma  whose 


wavalangtha  art  within  an  ordar  of  magnituda  of  tha  objact 
dimanaion  ia  damonatratad.  Tha  baaic  mathoda  diacuaaad 
hava  baan  pravioualy  idantifiad  in  tha  litaratura  aa 
natural  raaonanca  aatimation,  radar  imaging  from  ramp 
raaponaa  aignaturaa,  and  low  fraquancy  claaaification. 

Ralatad  common  faaturaa,  additional  inaight,  and  aoma  naw 
raaulta  ara  givan.  Applicationa  of  all  thraa  mathoda  to 
objacta  ranging  from  almpla  to  complax  gaomatriaa  ara 
daacrlbad. 

Murphy,  W.D.  and  M.L.  Bamaba,  "Numarical  Procaduraa  for  Solving 

Nonaymmatric  Eiganvalua  Problama  Aaaociatad  with  Optical  Raaonatora" 

Applied  Optica.  Vol.  17,  No.  IS,  pp.  2358-2365,  1  Auguat  1978. 

Tha  Prony  mathod  ia  axtandad  to  handla  tha  nonayrmnatric 
algabraic  aiganvalua  problam  and  improvad  to  aaarch  auto¬ 
matically  for  tha  numbar  of  dominant  aiganvaluaa.  A 
aingla  itarativa  algorithm  ia  givan  to  computa  tha 
aaaociatad  aiganvactora.  Raaolution  atudiaa  uaing  tha  QR 
mathod  ara  mada  in  ordar  to  datarmina  tha  accuracy  of  tha 
matric  approximation.  Numarical  raaulta  ara  givan  for  both 
almpla  wall-daflnad  raaonatora  and  mora  complax  advancad 
daaigna  containing  multiple  propagation  gaomatriaa  and 
mlaalignad  mirror a. 


McDonough,  R.N. ,  "Representation  and  Analysis  of  Signals,  Part  X V:  Matchsd 
Exponents  for  the  Rsprsssntation  of  Signals,"  Ths  Johns  Hopkins  Univsrsity, 
Department  of  Electrical  Engineering,  Baltimore  Maryland,  30  April  1963. 

Approximation  by  sums  of  exponentials  is  considered  from 
various  points  of  view,  such  as  the  time-domain  approximation 
problem  of  network  synthesis,  in  which  an  analytic  time 
function  is  to  be  approximated,  or  situations  in  which  dis¬ 
crete  measured  values  of  some  dynamic  process  are  to  be 
smoothed  and  fitted  using  an  exponential  model.  Prony's 
method  for  approximate  determination  of  the  exponents  is 
exhaustively  reviewed  in  all  its  forms. 

McDonough,  R.N.  and  W.H.  Huggins,  "Best  Least-Squares  Representation  of 
Signals  by  Exponentials,"  IEEE  Transactions  on  Automatic  Control.  Vol. 
AC-13,  Ho.  4,  pp.  408-412,  August  1968. 

In  this  paper  the  approximation  of  a  given  real  time  function 
over  (0,»)  by  a  linear  combination  of  a  given  number  n  of 
exponentials  is  considered,  such  that  the  Integrated  squared 
error  is  minimized  over  both  the  n  coefficients  of  the  linear 
combination  and  the  n  exponents  used.  The  usual  necessary 
condition  for  statlonarlty  of  the  integrated  squared  error 
leads  to  a  set  of  2n  simultaneous  equations,  nonlinear  in  the 
exponents.  This  condition  is  Interpreted  in  the  geometric 
language  of  abstract  vector  spaces,  and  an  equivalent  con¬ 
dition  Involving  only  the  exponents,  with  the  coefficients 
suppressed,  is  developed.  It  is  next  indicated  how  this 
latter  condition  can  be  applied  to  signals  which  are  not 
known  analytically,  but  only,  for  example,  as  voltages  re¬ 
corded  on  magnetic  tape,  or  as  a  table  of  sampled  values. 

The  condition  still  in  effect  requires  solution  of  nonlinear 
algebraic  equations,  and  a  linear  Iterative  method  is  pro¬ 
posed  for  this  purpose.  Finally,  the  procedure  Is  illustrated 
with  a  simple  example. 


Osborne,  M.R.,  "A  Class  of  Nonlinear  Regression  Problems",  in  Data 
Representation.  University  of  Queensland  Press,  pp.  94-100,  1970. 


This  paper  considers  the  problem  of  determining  the  para¬ 
meters  a,  g  in  the  model  equation 
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made  at  a  sequence  of  distinct  time  point  t. ,  t-  ...,  t  • 
It  is  assumed  that  the  e.  (the  experimental1errofs)  are  n 
independent,  normally  distributed,  random  variables  with 
mean  zero  and  standard  deviation  a  so  that  least  squares 
methods  are  appropriate.  It  is  also  convenient  to  assume 
that  the  functions  f  (t,  0)  form  a  Chebyshev  set  for  the 
values  of  g  of  interest.  “  This  means  that  no  nontrivial 
linear  comSinatlon  of  the  f  can  vanish  identically  on  any 
set  of  k  distinct  values  of1t.  It  is  convenient  to  refer 
to  a  regression  problem  of  the  above  form  as  separable,  and 
refer  to  a  and  g  as  the  linear  and  nonlinear  parameters 
respectively. 
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Parsons,  D.H. ,  "Biological  Problems  Involving  Sums  of  Exponential 
Functions  of  Time:  A  Mathematical  Analysis  that  Reduces  Experimental  Time" 
Mathematical  Biosciences .  ft 2,  pp.  123-128,  1968. 

The  analysis  given  in  this  article,  which  is  based  on  the 
theory  of  difference  equations,  enables  one  to  determine 
the  integer  rank  in  advance,  and  then  determine  the  2p 
unknown  constants  by  means  of  simultaneous  linear  equations. 

Paul,  A.K. ,  "Anharmonic  Frequency  Analysis",  Mathematics  of  Computation. 
Vol.  26,  No.  118,  pp.  437-447,  April  1972. 

A  new  numerical  method  of  frequency  analysis  is  described, 
designed  mainly  to  search  for  discrete  frequencies  in  a 
time  series.  An  integral  transform  is  applied  twice  to  the 
data  for  different  reference  times.  A  complex  amplitude 
within  a  selected  narrow  frequency  band  is  obtained  for  each 
transform.  The  frequency  is  then  determined  from  the  phase 
change  of  the  complex  amplitude  over  the  different  of  the 
two  reference  times.  Very  high  precision  is  obtained,  which 
is  demonstrated  in  two  examples. 

Pearson,  L.W. ,  M.L.  Van  Blaricum,  and  R.  Mittra,  "A  New  Method  for  Radar 
Target  Recognition  Based  on  the  Singularity  Expansion  for  the  Target", 
Proceedings  of  the  1975  IEEE  International  Radar  Conference.  Arlington, 
Virginia,  pp.  452-456,  21-23  April  1975. 

The  results  of  an  initial  study  of  the  use  of  the  recently 
introduced  Singularity  Expansion  Method  (SEM)  characterization 
of  a  scatterer  as  a  means  of  radar  target  identification  are 
reported.  Three  significant  advantages  of  this  method  over 
target  identification  schemes  are  explored.  These  are  1)  the 
aspect-dependence  of  transient  backscattered  waveforms  can  be 
suppressed,  2)  the  transmitted  waveform  requirements  are  only 
moderately  demanding  on  hardware,  and  3)  the  SEM  poles  consti¬ 
tute  a  smaller  collection  of  numbers  to  be  presented  to 
pattern  recognition  algorithms.  Examples  are  presented  indi¬ 
cating  the  following  features  of  the  method:  viability  in  the 
presence  of  noise;  excitability  of  a  reasonable  number  of 
natural  frequencies  with  a  moderate  bandwidth  pulse;  and  the 
ability  of  a  few  poles  to  identify  a  target  from  among  a  col¬ 
lection  of  known  targets. 

Pearson,  L.W.  and  0.  R.  Roberson,  "The  Extraction  of  the  Singularity 
Expansion  Description  of  a  Scatterer  from  Sampled  Transient  Surface 
Current  Response",  IEEE  Transactions  on  Antennas  and  Propagation. 

Vol.  AP-28,  No.  2,  pp.  182-190,  March  1980. 
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A  method  is  presented  whereby  one  cen  extract  the  singu¬ 
larity  expansion  method  (SEM)  description  of  an  object's 
electromagnetic  scattering  response  from  spatially  sampled 
transient  surface  currents.  The  currents  are  excited  by  a 
known  excitation.  The  SEM  data  are  recoverable  to  the 
degree  that  spatial  coupling  and  (frequency)  spectral 
intensity  excite  a  given  SEM  mode.  Results  of  a  numerical 
study  of  the  method  using  the  transient  response  of  a  thin 
wire  are  reported.  The  data  used  in  the  study  were 
computed  using  a  time  domain  Integral  equation  technique. 

The  ultimate  utility  of  the  method  lies  in  the  recovery  of 
SEM  data  from  measured  data,  thus  admitting  complex-shaped 
objects  into  the  realm  of  SEM  description.  The  method  is 
based  on  a  Prony-type  pole /residue  extraction  procedure. 

Pearson,  L.W.  and  D.R.  Roberson,  "Extrapolation  of  Induced  Transient 
Surface  Currents  on  a  Scatterer  to  General  Band-Limited  Excitation 
Conditions",  University  of  Kentucky  Electromagnetics  Research  Report, 
pp.  1-49,  October  1978. 

A  technique  for  the  extraction  of  an  approximate  singularity 
expansion  description  of  the  electromagnetic  response  of  a 
scatterer  is  presented.  The  description  is  extracted  from 
transient  waveforms  which  are  spatial  samples  of  the 
transient  current  on  the  scatterer  excited  by  a  known 
excitation.  The  method  can  use  any  reliable  Prony-type 
processor  for  exponential  analysis  of  the  transient  wave¬ 
forms  as  long  as  it  provides  physlcally-maaningful  poles 
from  the  waveforms.  The  fact  that  singularity  expansion 
poles  are  common  to  all  waveforms  on  the  object  provides  a 
computational  advantage  over  a  waveform-by-waveform  Prony- 
type  analysis.  A  "consensus  pole  set"  procedure  is 
described  to  exploit  this  redundancy.  Soma  Issues  relating 
to  the  completeness  of  a  set  of  transient  data  to  the  class 
of  excitations  to  which  the  singularity  expansion  is  applied 
is  discussed.  The  one-dimensional  example  of  a  thin -wire 
scatterer  is  used  to  test  the  practicability  of  the  method. 

The  ultimate  utility  of  the  method  is  demonstrated  by  way 
of  expansion  of  the  extracted  data  to  new  excitations. 
Recommendations  are  made  toward  an  implementation  of  the 
method  using  measured  data. 

Perkins,  J.F.  and  Charles  Cason,  "Effects  of  Small  Misalignments  in  Empty 
Unstable  Resonators",  Applied  Physics  letters.  Vol.  31,  No.  3,  pp.  198- 
200,  1  August  1977. 

Effects  of  misalignment  of  rectangular-cross-section 
unsteble  resonators  are  investigated  by  calculations  of 
eigenvalues  and  shapes  of  transverse  modes  of  empty  strip 
resonators.  Small  misalignments  can  substantially  Affect 


mode-loss  separation.  A  design  criterion  is  suggested  for 
.choice  of  Fresnel  numbers  such  that  effects  of  misalign¬ 
ments  are  minimized . 

Perl,  W. ,  "A  Method  for  Curve-Fitting  by  Exponential  Functions", 
International  Journal  of  Applied  Radiation  and  Isotropes ,  Vol.  8,  pp.  221- 
222. 

A  method  is  presented  for  fitting  exponential  functions  to 
data  of  the  decay  type.  In  an  elementary  "peeling-off" 
form  of  the  method,  the.  local  ordinate  is  plotted  against 
the  local  first  derivative  of  the  data.  The  straight  line 
portion  of  this  graph  at  large  time,  if  it  exists,  yields 
simultaneously  a  constant  (asymptotic  value)  plus  an 
exponential,  which  can  be  subtracted  from  the  given  curve 
and  the  procedure  repeated.  The  method  permits  judging  the 
data  in  simple  fashion  for  suitability  of  fit  by 
exponentials  simultaneously  with  the  fitting  itself.  The 
method  is  formulated  as  an  analytic  version  of  Prony's 
method,  which  permits  clarification  of  certain  difficulties 
of  the  latter.  Procedures  are  suggested  for  fitting 
exponentials  with  closely  spaced  half-lives,  and  for 
fitting  functions  other  than  exponential.  The  method  is 
illustrated  by  application  to  data  on  the  time  course  of 
respiratory  uptake  of  a  highly  fat-soluble  inert  gas  by 
human  subjects. 

Pogglo,  A.J.,  M.L.  Van  Blarlcum,  E.K.  Miller,  and  R.  Mlttra,  "Evaluation 
of  a  Processing  Technique  for  Transient  Data",  IEEE  Transactions  on 
Antennas  and  Propagation.  Vol.  AP-26,  No.  1,  pp.  165-173,  January  1978. 

A  data  processing  technique,  namely  Prony's  method,  is 
evaluated.  The  U3e  of  the  algorithm  for  EMP  and  other 
transient  problems  is  illustrated  and  the  relationship 
between  the  waveform  parameterization  and  the  singularity 
expansion  method  (SEM)  is  described.  The  impulse  response 
of  a  synthesizer  network  is  obtained  using  Prony  processing 
of  the  output  of  the  double  exponential  excited  circuit. 

Difficulties  which  arise  in  this  type  of  processing  such 
as  rank  deficiency,  aliasing,  and  noise  effects  are  con¬ 
sidered  and  methods  for  alleviation  such  as  filtering  and 
rank-overspecification  are  introduced  and  evaluated.  The 
overall  status  of  waveform  parameterization  as  applied  to 
SEM  is  considered  and  suggestions  for  future  research  pre¬ 
sented  . 

Prewitt,  J.F.  and  K.S.  Kunz,  "SEM  Pole  Extraction:  A  Feasibility  Study", 
Mission  Research  Corp.,  AMRC-N-61,  March  1978. 
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This  communication  establishes  typical  experimental 
electromagnetic  transient  response  measurement  limitations, 
and  applies  these  limitations  to  the  natural  mode 
decomposition  of  the  response  of  a  variety  of  scatterers 
due  to  EMP. 

Price,  H.J.,  "An  Improved  Prony  Algorithm  for  Exponential  Analysis",  AFWL 
Mathematics  Motes  #59,  November  1978. 

An  improved  Prony  algorithm  is  described  which  produces 
a  lower  sum  of  squares  of  residuals  than  the  usual  Prony 
algorithm.  Indeed,  the  new  algorithm  produces  the  lowest 
sun  of  squares  of  residuals  possible.  A  simple  numerical 
example  is  used  to  compare  the  new  and  old  algorithms. 

Within  the  context  of  the  new  algorithm,  a  technique  to 
constrain  poles  is  also  described. 

Price,  M.L.,  et  al.,  "Electrical  and  Photon  Tests  of  a  Resonant  Satellite 
Shape,"  IEEE  Transactions  on  Nuclear  Science.  Vol.  NS-25,  No.  6,  pp.  1358- 
1364,  December  1978. 

A  resonant  body  has  been  tested  and  analyzed  to  determine 
its  SGEMP  response.  Both  electrical  and  photon  excitation 
were  used  and  their  results  compare  favorably  with  pre¬ 
dictions.  A  modal  analysis  of  the  experimental  data  using 
Prony's  algorithm  yielded  the  dominant  frequencies,  damping 
and  excitation  amplitudes.  Within  the  accuracy  of  the 
measurements,  the  presence  of  the  phot on- induced  electron 
cloud  had  no  effect  on  the  period  or  damping  times  of  the 
RESMOD  modes.  The  proximity  of  an  electrical  pulsar  and 
its  200-0hm  terminating  resistor  markedly  decreased  the 
damping  time  during  electrical  excitation. 

Prony,  R. ,  "Essai  Experimental  et  Analytique  sur  les  Lois  de  la  Dilatabilite 
de  Fluides  Elastlques  et  sur  Celles  de  la  Force  Expansive  de  la  Vaperu  de 
l'Alkool,  a  Differentes  Temperatures,"  J.  l'Scole  Polvtech.  (Paris), 

Vol.  1,  No.  2,  pp.  24-76,  1795. 

This  paper  has  been  abstracted  in  the  Section  2.0  of  this 
bibliography.  This  is  Prony's  original  paper. 


Quasi,  A.H.  and  F.R.  Spitznogle,  "Time-Limited  Signal  Representation  and 
Analysis  by  Using  Non-Fourier  Techniques,"  U.S.  Navy  Journal  of  Under¬ 
water  Acoustics,  Vol.  18,  No.  3,  pp.  397-402,  July  1968. 

Time-limited  signals  are  represented  and  analyzed  by  using 
non-Fourier  techniques.  Results  show  that  non-Fourier 
techniques  can  be  superior  to  the  conventional  Fourier 
approach  for  the  frequency  analysis  of  time  limited  signals. 


Salzer,  H.E. ,  "Some  Extensions  of  Prony  Approximation",  Zeitschrift  fur 
Angewandte  Mathematik  und  Mechanik.  57,  #4,  pp.  269-271,  1977. 


Prony1 e  method  of  approximating  a  function  f(s)  by  a  sum 

n  a^x 

of  exponential  terms  A^  a  1  ,  where  both  a^  and  A^ 

are  to  be  determined,  is  extended  to  oscillatory  interpola¬ 
tion  of  any  order,  and  also  to  the  direct  global  approxima¬ 
tion  to  the  solution  of  the  linear  differential  equation 

eJL o  0r(x)  y(r)<x) "  f(*>  *°* 

Sarkar,  T.K. ,  J.  Nebat,  and  D.D.  Weiner,  "Suboptimal  System  Approximation/ 
Identification  with  Known  Error",  Air  Force  Weapons  Laboratory,  AFWL-TR- 
77-200,  September  1977.  Also  Mathematics  Motes  <1149. 

This  paper  presents  a  noniterative  method  for  approximating 
empirical  signals  over  [0,*]  by  a  linear  combination  of 
exponentials.  The  technique  results  in  a  suboptimal 
approximation.  Motably,  the  dependence  of  the  suboptimal 
exponents  s^  on  the  integral  square  error  c  is  such  that 

lim  (e-0  s^  “  s^,  the  optimal  exponents.  The  method  may 

also  be  used  for  system  identification.  It  la  especially 
useful  when  the  system  is  modelled  by  a  black  box  and  one 
has  access  only  to  the  input  and  output  terminals  of  the 
system.  A  technique  is  demonstrated  to  find  the  multiple 
poles  of  a  system  along  with  the  residues  at  the  poles  when 
the  output  of  the  system  to  a  known  input  is  given.  Among 
the  advantages  of  the  method  are  its  natural  insensitivity 
to  noise  in  the  data  and  the  explicit  determination  of  the 
signal  order.  Representative  computations  are  made  of  the 
poles  from  the  transient  response  of  a  conducting  pipe 
tested  at  the  ATHAMAS-I  EMP  simulator. 

Sarkar,  T.K. ,  V.K.  Jain,  J.  Mebat,  and  D.D.  Weiner,  "A  Comparison  of  the 
Pencil-o f-FunctJ on  Method  with  Prony 's  Method,  Weiner  Filters  and  Other 
Identification  Techniques",  Rochester  Institute  of  Technology,  AFWL 
Mathematics  Motes  # 54,  December  1977. 

The  problem  of  interest  is  to  identify  the  transfer 
function  of  a  system  by  its  poles  and  residues  when  the 
noise  contaminated  input  and  output  are  specified.  The 
first  aim  of  this  paper  is  to  illustrate  that  several 
different  formulations  for  characterizing  the  impulse 
response  of  a  system  yield  the  same  set  of  poles  as  is 
obtained  in  the  case  of  a  discrete  Weiner  filter.  The 
second  goal  is  to  show  how  different  formulations 
regularize  the  ill-posed  system  identification  problem. 


Ic  is  demonstrated  that  the  Weiner  filter  la  not  always 
realizable  aa  a  causal  rational  function.  When  the  order 
of  the  filter  la  specified  a  priori,  the  resulting  filter 
may  no  longer  be  optimum.  Finally,  representative 
computations  are  made  of  the  poles  from  the  transient 
response  of  a  conducting  pipe  tested  at  the  ATHAMAS-I  EMP 
simulator  to  demonstrate  the  stability,  reliability, 
consistency,  and  accuracy  of  the  results  obtained  by  the 
pencil-of-f unction  method. 

Sarkar,  T.K. ,  D.D.  Weiner,  and  V.K.  Jain,  "Some  Mathematical  Considerations 
in  Dealing  with  the  Inverse  Problem",  IEEE  Transactions  on  Antennas  and 
Propagation.  Vol.  AP-29,  No.  2,  pp.  373-379,  March  1981. 
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Many  problems  of  mathematical  physics  can  be  formulated  in 
terms  of  the  operator  equation  Ax  ■  y,  where  A  is  an 
lntegro-differential  operator.  Given  A  and  x,  the  solution 
for  y  is  usually  straightforward.  However,  the  inverse 
problem  which  consists  of  the  solution  for  x  when  given  A 
and  y  is  much  more  difficult.  The  following  questions 
relative  to  the  Inverse  problem  are  explored.  1)  Does 
specification  of  the  operator  A  determine  the  set  {y}  for 
which  a  solution  x  is  possible?  2)  Does  the  Inverse  problem 
always  have  a  unique  solution?  3)  Do  small  perturbations 
of  the  forcing  function  y  always  result  in  small  perturba¬ 
tions  of  the  solution?  4)  What  are  some  of  the  considera¬ 
tions  that  enter  into  the  choice  of  a  solution  technique 
for  a  specific  problem?  The  concept  of  an  ill-posed 
problem  versus  that  of  a  well-posed  problem  is  discussed. 
Specifically,  the  manner  by  which  an  ill-posed  problem  may 
be  regularized  to  a  well-posed  problem  is  presented.  The 
concepts  are  illustrated  by  several  examples. 

Scharf,  L.L.  and  A. A.  Beex,  "Covariance  Sequence  Approximation  for 

Parametric  Spectrum  Modelling",  Technical  Report  August  80-1,  Colorado 

State  University,  Ft.  Collins,  Colorado  80523. 


Parametric  methods  of  spectrum  analysis  are  founded  on 
finite-dimensional  models  for  convariance  sequences. 
Rational  spectrum  approximants  for  continuous  spectra  are 
based  on  autoregressive  (AR) ,  moving  aversge  (MA)  or  auto¬ 
regressive  moving  average  (ARMA)  models  for  covariance 
sequences.  Line  spectrum  approximants  to  discrete  spectra 
are  based  on  cosinusoidal  models  for  covariance  sequences. 
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In  this  paper  the  authors  make  the  point  that  a  wide  CJ 

variety  of  spectrum  types  admit  to  modal  analysis  wherein 

the  modes  are  characterized  by  amplitudes,  frequencies, 

and  damping  factors.  The  associated  modal  decomposition 

is  appropriate  for  both  continuous  and  discrete  components 
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of  the  spectrum.  ...a  domain  of  attraction  for  tha 
dacompoaition  includaa  ARMA  sequences,  harmonically-  or 
nonharmonically-rolated  sinusoids,  danspad  sinusoids,  whita 
noise,  and  linear  combinations  of  these. 

The  parametric  spectrum  analysis  problem  now  becomes  one  of 
identifying  mode  parameters.  This  is  achieved  by  solving 
two  modified  least-squares  problems.  Numerical  results  are 
presented  to  illustrate  the  identification  of  mode 
parameters  and  corresponding  spectra  from  finite  records  of 
perfect  and  aatimatad  covariance  sequences.  The  results 
for  sinusoids  and  sinusoids  in  white  noise  are  interpreted 
in  terms  of  inphase  and  quadrature  effects  attributable  to 
the  finite  record  length. 

Scharf,  L.L.,  A. A.  Beex  and  T.  von  Reyn,  "Modal  Decomposition  of 
Covariance  Sequences  for  Parametric  Spectrum  Analysis",  contributed  paper, 
1981  International  Conferenca  on  Acoustics,  Speech,  and  Signal  Processing, 
Atlanta,  March  30,  1981. 

In  this  paper  the  authroa  make  the  point  that  a  wide 
variety  of  spectrum  types  admit  to  modal  analysis  wherein 
the  modes  are  characterized  by  amplitudes,  frequencies, 
and  damping  factors.  The  associated  modal  decomposition  is 
appropriate  for  both  continuous  and  discrete  components  of 
the  spectrum.  The  domain  of  attraction  for  tha  decomposi¬ 
tion  includes  ARMA  sequences,  harmonically-  or  nonharmoni- 
cally-relatad  sinusoids,  damped  sinusoids,  white  noise,  end 
linear  combinations  of  these. 

Numerical  results  are  presented  to  illustrate  the  identifi¬ 
cation  of  mode  parameters  and  corresponding  spectra  from 
finite  records  of  perfect  and  estimated  covariance  sequences. 

The  results  for  sinusoids  and  sinusoids  in  white  noise  are 
interpreted  in  terms  of  lnphase  and  quadrature  affects 
attributable  to  the  finite  record  length. 

Schaubert,  D.H. ,  "Application  of  Prony's  Method  to  Time-Domain 
Ref lactometer  Data  and  Equivalent  Circuit  Synthesis",  HDL-TR-1857,  Harry 
Diamond  Laboratories,  June  1978.  Also  IEEE  Transactions  on  Antennas  and 
Propagation.  Vol.  AP-27,  No.  2,  pp.  180-184,  March  1979. 

The  singularity  expansion  of  time-domain  reflectomater  data 
has  been  computed  and  used  to  calculate  lumped-element 
equivalent  circuits  for  the  impedance  of  some  typical 
antennas.  Prony's  algorithm  was  used  to  obtain  the  poles 
of  the  antennas'  terminal  voltage  waveform  due  to  a  step- 
like  excitation.  The  residues,  however,  were  calculated 
subject  to  the  unlform-srror  uorm  instead  of  the  least- 
squerad-error  norm.  Although  it  requires  much  longer  to 
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calculate,  the  uniform-arror  norm  approximation  la  often 
orders  of  magnitude  better  than  the  least-squared-error 
approximation.  The  physical  realizability  of  impedance 
functions  obtained  from  experimental  data  has  been 
investigated.  Also,  the  procedures  necessary  to  synthesize 
the  lumped-element  network  have  been  evaluated.  A  simple 
partial-fraction  expansion  of  a  realizable  impedance 
function  representing  an  antenna  was  found  to  yield 
nonrealizable  element  values.  A  general  network  synthesis 
procedure,  such  as  Bruns' s  method,  is  required. 

Schultz,  B.R. ,  "Estimation  of  Pulse  Transfer  Function  Parameters  by  Qua. 

linearization",  IEEE  Transactions  on  Automatic  Control,  pp.  424-426, 

August  196S. 

An  iterative  algorithm  for  the  identification  of  single- 
input  single-output  linear  stationary  discrete  system  is 
developed  using  the  method  of  quaailinearlzatlon.  The 
resulting  procedure  is  similar  to  mode  1  of  the  method  of 
Stelglitz  and  McBride  but  has  the  advantage  of  the 
quadratic  convergence  property  of  quasllinearlzatlon.  It 
is  shown  that  this  algorithm  becomes  mode  1  if  the 
measured  plant  output  is  used  in  the  calculations  in  place 
of  the  model  output.  Consequently,  the  two  methods  are 
extremely  compatible  and  it  is  a  simple  matter  to  combine 
them  in  a  single  program,  which  generates  its  own  initial 
estimate,  has  the  wide  range  of  convergence  of  mode  1,  and 
possesses  the  quadratic  convergence  property  of  quasi¬ 
linearization  for  final  convergence  to  a  solution.  The 
method  also  permits  the  estimation  of  plant  initial 
conditions  in  those  cases  whare  they  must  be  considered. 

Results  of  a  few  numerical  applications  are  discussed. 

Scrivner,  G. J. ,  "Prony's  Analysis  in  the  Presence  of  Noise",  Computer 

Sciences  Corporation,  Report  II/3-C-235,  October  1977.  Also  AFWL 

Mathematical  Note  // 67. 

This  report  studies  a  methodology,  based  on  Prony's 
algorithm,  for  extracting  complex  frequencies  and 
associated  residues  directly  from  transient  response  data. 

The  methodology  is  not  new  but  has  been  receiving 
increasing  attention  due  to  its  potential  usefulness  in  a 
variety  of  different  applications  in  quits  diverse  scien¬ 
tific  disciplines.  Although  this  particular  study  was  to 
evaluate  the  potential  usefulness  of  Prony's  method  for 
analysis  of  System  Generated  Electromagnetic  Pulse  (SGEMP) 
experimental  data,  the  results  obtained  are  of  a  much 
broader  nature.  The  report  presents  a  classical  derivation 
of  the  matrix  equation  which  result  if  one  uses  Prony's 
algorithm  in  conjunction  with  a  least-squares  criterion. 


The  resulting  matrix  is  studied  with  respect  to  tho  nature 
of  its  eigenvalue  structure  and  an  error  analysis  is 
developed  which  highlights  the  Importance  of  this  structure. 

The  influence  of  noise  on  the  eigenvalue  structure  and  its 
Impact  on  the  determination  of  the  number  of  polee  in  the 
data  are  examined.  Optimisation  concepts  relative  to 
advantageous  modification  of  ths  eigenvalue  structure  of 
the  Prony  matrix  !*:e  qualitatively  discuased.  Difficulties 
associated  with  the  two  state  applicetlon  of  the  least- 
squares  method  in  conventional  Prony  analysis  are  cited  and 
an  iterative  method  of  removing  this  shortcoming  described. 

Numerical  calculations  illustrating  the  majority  of  the 
above  concepts  are  prssented  and  discussed. 

S legman,  E.  and  H.Y.  Miller,  "Unstable  Optical  Resonator  Loss  Calculations 
,  Using  the  Prony  Method",  Applied  Optica.  Vol.  9,  No.  12,  pp.  2729-2736, 

December  1970. 

The  eigenvalues  for  all  the  significant  low-order  resonant 
modes  of  an  unstable  optical  resonator  with  circular 
mirrors  are  computed  using  an  eigenvalue  method  called  the 
t  Prony  method.  A  general  aqui valence  relation  is  alao 

given,  by  means  of  which  one  can  obtain  the  design 
paramsters  for  a  single-ended  unstable  resonator  of  the 
type  usually  employed  in  practical  lasers,  from  the 
calculated  or  tabulated  values  for  an  equivalent  symmetric 
or  double-ended  unstable  resonator. 

\ 

Simon,  W. ,  "A  Method  of  Exponential  Separation  Applicable  to  Small 
Computers",  Phvs.  Mad.  Biol..  Vol.  15,  No.  2,  pp.  353-360,  1970. 

A  method  has  bean  developed  for  separation  of  the  components 
of  multlcompertment  exponentiels  which  operates  on  a  small 
i  computer.  It  seams  free  of  many  of  the  difficulties  of 

other  methods  and  provides  e  reasonable  test  for  the  number 
of  components. 

Smith,  W.D.  and  D.L.  Lager,  "Parametric  Characterization  of  Random 
Processes  Using  Prony* 8  Method",  Lawrence  Livermore  Laboratory,  UCRL- 

>  52673,  25  May  1979. 

The  compact  parametric  characterization  of  a  random 
process  is  often  valuable  in  problems  Involving  signal 
classification  and  system  identification.  The  method  of 
Prony  suggest  two  different  approaches  to  obtaining  such  a 

>  characterization  from  a  finite  record  of  sampled  date. 

The  first  approach  is  to  extract  the  parameters  of  ths 
estimated  autocorrelation  waveform.  The  second  approach  is 
to  estimate  the  parameters  of  a  white  noise  driven  linear 
system  that  generates  an  output  with  the  same  character¬ 
istics  as  the  given  random  process. 


These  cwo  approaches  ars  compared  with  other  common  signal 
processing  methods.  Also,  the  performance  of  each  approach 
is  evaluated  theoretically  and  with  computer-simulated  and 
experimentally  recorded  data.  Each  approach  Is  found  to 
rapidly  characterize  both  s  ationary  and  nonstatlonary  random 
processes.  The  first  approach  is  superior  for  the  accurate 
paramatric  characterization  of  a  stationary  process.  Alter¬ 
nately,  the  more  precise  second  approach  is  beat  for  monitoring 
changes  with  time  in  random  process  parameters. 

Soussou,  J.E.,  F.  Moavenzadeh,  and  M.H.  Gradowczyk,  "Application  of  Prony 
Series  to  Linear  Vaicoelaaticity",  Transactions  of  the  Society  of 
Rheology.  Vol.  14,  No.  4,  pp.  373-584,  1970. 

A  new  method  for  description  of  viscoelastic  functions  with 
Prony  series  of  exponentials  is  presented  and  compared  with 
existing  methods.  This  method  places  constraints  on  the 
coefficients  of  the  series  which  guarantees  smoothness  of 
thu  functions  and  a  discrete  spectral  representation.  An 
optimization  technique  is  used  to  obtain  the  solution.  The 
resulting  series  permit  the  use  of  Whittaker's  method  of 
solution  of  Integral  equations  which  is  specialized  to  the 
use  of  viscoelastic  analysis,  and  an  exact  relationship  is 
obtained  between  creep  and  relaxation  functions. 

Spitznogle,  F.R.  and  A.H.  Quazl,  "Representation  and  Analysis  of  Tiwa- 
Llmitad  Signals  Using  a  Complex  Exponential  Algorithm",  Journal  of  the 
Acoustical  Society  of  America.  Vol.  47,  No.  5,  pp.  1150-1155,  May  1970. 

A  complex  exponential  algorithm  developed  for  the 
representation  and  analysis  of  time-limited  signals  is 
defined,  and  its  evaluation  with  respect  to  conventional 
discrete  Fourier  techniques  is  discussed.  It  is  shown  that 
for  a  given  length  of  a  signal  containing  discrete  frequency 
information,  sampled  at  least  to  the  Nyquist  criterion,  the 
complex  exponential  algorithm  can  often  provide  increased 
frequency  resolution  over  standard  Fourier  techniques.  It 
is  also  shown  that  the  complex  exponential  algorithm 
provides  an  improved  mechanism  over  Fourier  techniques  for 
interpolation  between  points  in  a  sampled  signal  containing 
discrete  frequency  components  in  the  presence  of  board-band 
noise.  The  effects  of  noise  in  the  complex  exponential 
technique  and  the  computational  difficulties  associated  with 
the  present  complex  exponential  algorithm  are  also  discussed. 


Texas  Instruments  Incorporated,  "Annual  Report  -  Representation  and 
Analysis  of  Sonar  Signals  Volume  I:  Improvements  in  the  Complex  Ex¬ 
ponential  Signal  Analysis  Computational  Algorithm,"  Texas  Instruments, 
Inc.,  Dallas,  Tf -.as,  31  March  1970. 

This  volume  of  the  1969-70  Annual  Report  on  Contract 
N00014-69-C-0315  describes  the  results  of  a  study  to  re¬ 
duce  the  computer  memory  and  computational  time  required 
to  obtain  a  complex  exponential  signal  representation. 

As  described  in  the  report,  substantial  reductions  in  com¬ 
putational  complexity  have  been  achieved  through  improvements 
in  the  Prony  method  for  solving  the  nonlinear  system  of 
equations  associated  with  a  complex  exponential  signal 
analysis.  Projections  of  the  real-time  capability  of  a 
hard-wired  complex  exponential  processor  using  state-of- 
the-art  hardware  are  given.  In  addition,  two  iterative 
complex  exponential  algorithms  based  on  the  Newton-Raphson 
and  gradient  techniques  for  solving  a  nonlinear  system  of 
equations  are  described. 

Trivett,  David  H.,  "A  Modified  Prony  Method  Approach  to  Echo-Reduction 
Measurements  of  Time-Limited  Transient  Signals,"  NRL  Memorandum  Report 
4172,  June  16,  1980,  ADA  087211. 

A  modified  Prony  method  is  presented  for  measuring  the 
steady-state  echo  reduction  of  acoustic  panels.  The 
method  extrapolates  the  stead-state  amplitudes  from  the 
transient  portion  of  the  signal  allowing  time-limited 
measurements.  The  method  is  applied  to  measurements  of 
square  panels  76  cm  on  an  edge  and  0.95  cm- thick  steel 
and  aluminum  in  the  frequency  range  of  3-10  kHz.  The 
signals  were  time-limited  to  200ys  (0.6  -  2. OX)  by  the 
arrival  of  the  diffracted  signal  from  the  panel  edges. 

Results  are  compared  with  theoretical  values  and  a  pro¬ 
gram  listing  is  included. 


Van  Blarlcum,  M.L.  and  R.  Mittra,  "A  Technique  for  Extracting  the  Poles 
and  Roaiduea  of  a  System  Directly  from  its  Transient  Response/'  IEEE 
Transactions  on  Antennas  and  Propagation.  Vol.  AP-23,  No.  6,  pp.  777- 
781,  November  1975. 


A  novel  approach  for  systematically  deriving  the  complex 
poles  and  residues  of  a  structure  from  a  set  of  time-domain 
data  is  presented.  The  method  is  based  on  Prony's  algorithm 
which  involves  the  inversi.u  of  two  matrices  and  a  solution 
of  the  zeros  of  an  nth  degree  polynomial,  n  being  the  number 
of  desired  poles.  Two  numerical  examples  are  presented  and 
several  numerical  advantages  of  this  method  are  discussed. 

Van  Blarlcum,  M.L.,  "Techniques  for  Extracting  the  Complex  Resonances  of 
a  System  Directly  from  its  Transient  Response,"  Ph.D.  Dissertation,  De¬ 
partment  of  Electrical  Engineering,  University  of  Illinois,  December  1975; 
also  AFWL  EMP-Interaction  Note  3d. 

t 


The  singularity  expansion  method  is  a  techniques  used  for 
writing  the  transient  response  of  a  structure  as  a  sum  of 
exponentially  damped  sinusoids.  In  order  to  apply  this 
method,  it  is  first  necessary  to  obtain  the  singularities 
of  the  system  being  studied.  The  conventional  approach  for 
determining  the  singularities  of  a  system  is  based  on  an 
iterative  search  procedure  that  seeks  the  zeros  of  the 
system  determinant  in  the  complex  frequency  plane.  The 
alternative  approach  of  extracting  the  system  singularities 
directly  f  :om  the  transient  response  function  is  discussed. 
The  method  developed  is  baaed  on  Prony's  algorithm.  Some 
basic  problems  which  are  associated  with  the  use  of  Prony's 
algorithm  are  discussed  and  solutions  are  obtained. 


It  is  demonstrated  that  Prony's  method  is  applicable  to 
systems  with  multiple  as  well  as  simple  pole  singularities. 
Two  techniques  are  presented  for  systematically,  determining 
the  number  of  poles  contained  in  a  transient  response. 


Van  Blarlcum,  M.L.  and  R.  Mittra,  "PSPL:  A  Computer  Program  for 
Extracting  the  Poles  and  Residues  from  a  Transient  Signal  Using  Prony’s 
Method",  Electromagnetics  Laboratory  Scientific  Report  No.  7602, 
University  of  Illinois,  Urbana,  Illinois,  March  1976. 

This  report  describes  the  Fortran  program  PSPL,  which  Is 
used  to  extract  the  poles  and  residues  of  a  system  from  a 
set  of  discrete  transient  data  using  a  least-squares 
version  of  Prony's  method.  The  code  has  been  written  so 
that  multiple  as  well  as  simple  poles  can  be  obtained  if 
they  are  present  in  the  data.  The  listing  of  the  program 
given  In  this  report  Is  for  use  on  the  CDC-7600  computer 
but  can  be  easily  adapted  to  other  machines. 
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Van  Blarlcum,  M.L.  and  R.  Mittra,  "PEIGEN:  A  Computer  Program  for  the 
Determination  of  the  Number  of  Poles  in  a  Transient  Signal", 

Electromagnetics  Laboratory  Report  No.  76-5,  University  of  Illinois, 

Urbana,  Illinois,  June  1976. 

This  is  a  users  program  for  the  PEIGEN  program  which  is 
used  to  determine  the  number  of  poles,  N  in  a  set  of 
transient  data  assuming  that  the  data  are  made  up  of  a 
finite  sum  of  exponentials.  PEIGEN  is  written  as  a 
companion  to  the  program  PSPL,  which  is  used  to  extract  the 
N  poles  from  the  transient  data  once  the  value  of  N  has 
been  determined.  The  method  uses  the  eigenvalue  method. 

Van  Blarlcum,  M.L. ,  "The  Use  of  Prony's  Method  as  an  Analysis  Tool  Applied 
to  Electrical  Excitation  of  FltSatCom",  Mission  Research  Corporation  Report, 
MRC-R-3C1,  January  1977. 

This  report  discusses  the  use  of  the  numerical  algorithm 
known  as  Prony's  method  as  a  tool  for  the  analysis  of 
electrical  transient  experimental  data.  The  data  to  which 
these  techniques  are  applied  as  examples  are  results  of 
electrical  excitation  tests  of  a  structural  mockup  of  the 
FltSatCom  satellite. 

Van  Blarlcum,  M.L.  and  R.  Mittra,  "Problems  and  Solutions  Associated  with 
Prony's  Method  for  Processing  Transient  Data",  IEEE  Transactions  on 
Antennas  and  Propagation.  Vol.  AP-26,  No.  1,  pp.  174-182,  January  1978. 

Three  difficulties  associated  with  Prony's  method  are 
studied.  These  are  the  extension  of  the  method  to  allow 
multiple  poles,  the  development  of  techniques  for 
determining  the  number  of  poles  contained  in  the  transient 
data,  and  the  effects  of  noise  In  the  data  on  the  numerical 
procedures.  Solutions  to  these  difficulties  are  studied, 
and  numerical  samples  are  presented. 
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Weaver,  J.J.,  D.  Mg,  and  D.  L.  Lager,  "Analytical  Qualification  of  System 
Identification  (Modal  Analysis)  Codes  for  Use  in  the  Dynamic  Testing  of 
Nuclear  Power  Plant  Structure",  Lawrence  Livermore  Laboratory,  UCID- 18144, 
2  January  1980. 

This  report  describes  the  analytical  evaluation  of  two 
particular  system  identification  codes  used  at  Lawrence 
Livermore  Laboratory.  Both  codes  are  eigenparameter  iden¬ 
tification  codes,  however,  one  uses  a  time  domain  approach 
while  the  other  a  frequency  domain  approach.  The  evaluation 
was  accomplished  by  analytically  generating  several  time 
history  signals  in  which  the  true  modal  parameters  were  known. 

These  time  histories  ranged  from  widely  spaced  modes  with 
spacing  factors  of  6  percent.  These  signals  were  then 
polluted  with  various  levels  of  simulated  measurement  noise 
and  the  ability  of  our  computer  codes  to  extract  the  para¬ 
meters  from  this  noisy  data  was  evaluated. 

Weiss,  L.  and  R.N.  McDonough,  "Prony's  Method,  Z-Trans forms ,  and  Pade 
Approximation",  SIAM  Review.  Vol.  5,  No.  2,  pp.  145-149,  April  1963. 

The  purpose  in  this  note  is  to  point  out  that  the  original 
form  of  Prony's  method  may  also  be  regarded  as  a  Pade 
approximation  in  the  Z-transform  domain. 

Weygandt,  C.N.  and  N.N.  Puri,  "Identification  of  a  Linear  System  From 
Discrete  Values  of  Input-Output  Data",  1966  Proceeding  4th  Annual 
Allerton  Conference  on  Circuit  and  System  Theory,  po.  903-912,  1966. 

In  this  paper  the  authors  describe  a  method  for 
determining  the  transfer  function  of  a  physical  system 
having  period  input  and  output,  both  of  which  can  be 
measured.  The  method  depends  upon  taking  2n  equidistant 
samples  of  the  input  and  the  output  within  one  period,  n 
being  the  number  of  poles  the  system  transfer  function  is 
assumed  to  have.  These  samples  are  used  to  compute  (n+1) 
determinants  and  thereby  the  coefficients  and  the  roots  of 
an  nth  order  polynomial.  The  roots  of  this  nth  order 
polynomial  are  further  processed  to  obtain  the  poles  of  the 
transfer  function.  Once  the  location  of  poles  is 
determined,  a  very  simple  computer  routine  determines  the 
residues  of  these  poles  and  thereby  the  zeros  and  the  gain 
constant  are  determined.  An  essential  feature  of  the  scheme 
is  that  it  determines  the  poles  independent  of  the  zeros. 

In  practical  examples  the  number  of  samples  taken  IS  much 
larger  than  2n,  thereby  generating  a  large  number  of  nth 
order  polynomials.  The  average  of  the  roots  of  these 
polynomials  are  used  to  locate  the  poles. 
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White,  O.R.,  "Inversion  of  the  Limb-Darkening  Equation  Using  the  Prony 
Algorithm",  The  Aatophysical  Journal.  Vol.  152,  pp.  217-223,  April  1968. 

King's  interesting  application  of  the  Prony  algorithm  to 
the  inversion  of  the  limb-darkening  equation  is  modified 
and  generalized  so  that  it  can  be  applied  rigorously  to 
limb-darkening  curves,  spectral  line  profiles,  and 
multiplet  line  intensities.  Since  King's  formulation 
requires  physically  inaccessible  observations,  a  change  of 
variable  is  Introduced  in  the  limb-darkening  integral  to 
avoid  this  difficulty.  Calculations  on  noisy  data 
confirm  the  self-limiting  property  of  this  inversion 
method.  An  approximate  method  for  constructing  the  smooth 
S(t)  distribution  from  the  slab  solution  is  presented. 

Wiscombe,  W.J.  and  J.  W.  Evans,  "Exponential-Sum  Fitting  of  Radiative 
Transmission  Functions",  Journal  of  Computational  Physics.  24,  pp.  416- 
444,  1977. 

Fitting  transmission  functions  with  exponential  sums  is  the 
basis  for  a  widely  used  approximation  for  calculating 
spectrally  integrated  radiative  fluxes  in  planetary 
atmospheres,  expecially  when  both  line  absorption  and 
scattering  are  important.  The  error  in  this  method 
depends  crucially  on  the  accuracy  of  the  fits,  but 
unfortunately  exponential-sum  fitting  is  a  classical  ill- 
conditioned  problem  of  numerical  analysis.  Previous 
techniques  devised  for  exponential- sum  fitting  are  often 
unsatisfactory  in  this  application.  This  part  presents  a 
new  method  which  sidesteps  the  ill-conditioning,  guerentees 
convergence  to  the  unique  best  least-squares  fit,  gives 
positive  coefficients,  and  produces  fits  orders  of 
magnitude  more  accurate  than  any  which  have  so  far  been 
published.  The  method  is  demonstrated  to  be  capable  of 
recovering  an  exponential  sum,  given  data  sampled  from  that 
sum  and  rounded  to  as  few  as  two  decimal  places.  Sample 
fits  are  given  for  the  Goody  and  Malkmus  random  band 
models  and  the  Yamamoto  H2O  solar  absorption  data  in  order 
to  illustrate  the  high  accuracy  of  the  method.  The  affect 
on  the  fits  or  errors  in  the  transmission  data  is  examined 
in  some  depth. 

Wolf,  J.K. ,  "Decoding  of  Bose-Chaudhurl-Hocquenghem  Codes  and  Prony's 
Method  for  Curve  Fitting",  IEEE  Transactions  on  Information  Theory,  p.  608, 
October  1967. 

The  literature  is  filled  with  examples  of  a  common  set  of 
equations  which  arise  in  two  or  more  diverse  applications. 

The  purpose  of  this  correspondence  is  to  point  out  that 
such  a  situation  has  occurred  in  the  two  fields  of 


1)  algebraic  coding  theory,  and  2)  curve  fitting.  Two 
aspects  are  of  particular  Interest:  the  large  time  span 
which  separates  the  fundamental  works  in  these  two  areas, 
and  that  the  methods  of  solution  are  identical. 

Wynn,  P.,  "A  Note  on  the  Fitting  of  Certain  Types  of  Experimental  Data", 
Statistics  Neerlandica,  Vol.  16,  No.  2,  pp,  143-150,  1962. 

This  paper  is  concerned  with  the  problem  of  determining 
those  values  of  the  parameters  Z,  At,  ai(i-l,2, . . . ,h) 
which  provide  a  best  fit  of  the  form 

k  -a ,  t 

f*(t)  -  Z  +  L  A.e  1 
i-1 

to  the  observed  values  fg  -  f(tg) (s-1 ,2, . . . ,n)  in  the  sense 
k  9 

that  the  3um  Z  {fc‘*(t3)  -  f(tg)}  is  a  minimum. 
i-1 

Two  iterative  numerical  techniques  for  the  solution  of  this 
problem  are  discussed.  The  first,  a  variable  gradient  method, 
converges  rapidly  but  involves  a  relatively  large  amount  of 
computation;  the  second  involves  less  computation  but  con¬ 
verges  less  rapidly.  It  is  shown  that  by  application  of  a 
suitable  acceleration  technique  to  the  second  method,  its 
coo  'ergence  is  made  as  rapid  as  that  of  the  first. 
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Yengsc,  W.C.,  "Approximation  to  a  Specified  Time  Response,"  IRE  Trans¬ 
actions  on  Circuit  Theory,  pp.  152-163,  June  1962. 

This  paper  presents  a  procedure  by  which  specified  data  or  i 
function  of  time  h*(t)  can  be  approximated  by  trigonometric 
and/or  exponential  functions  of  time  h(t)  for  which  the 
Laplace  transformations  H(s)  are  known  and  can  be  expressed 
in  rational  fraction  form.  The  procedure  is  based  on  fitting 
h((t)  by  an  mth-order  difference  equation  whose  coefficients 
are  determined  by  a  least-squares  technique.  These  co¬ 
efficients  are  used  directly  to  determine  the  poles  of  H(s). 

The  zeros  of  H(s)  are  established  by  using  the  prescribed 
data  or  function  h*(t)  and  the  initial  value  theorem.  The 
approximate  function  of  time  is  obtained  by  taking  the  in¬ 
verse  Laplace  transformation  of  H(s) .  By  this  procedure 
not  only  is  an  approximation  obtained  for  h*(t)  in  the  time 
domain,  but  its  transform  is  also  found  in  rational  fraction 
form  suitable  for  realization  as  a  driving  point  or  transfer 
function.  Furthermore,  the  least-squares  technique  used  in 
determining  most  or  all  of  the  unknown  parameters  in  this 
procedure  tends  to  minimize  the  effect  of  random  errors  or 
noise  present  in  the  specified  data. 

Youla,  D.C.,  "The  Identification  of  Linear  Dynamical  Systems  from  Time- 
Domain  Measurements,"  Department  of  Electrical  Engineering  -  Farmlngdale, 
New  York  11735,  pp.  1-26,  Appendix  A1-A2,  1  September  1978. 

The  problem  of  identifying  the  Impulse  response  of  a  lumped, 
linear  time- invariant  system  from  a  discrete  set  of  its 
samples,  is  essentially  finite-dimensional  in  character. 

Prony's  method  is  distinguished  by  the  fact  that  it  re¬ 
cognizes  this  finite-dimensionality  from  the  outset  and 
thereby  succeeds  in  avoiding  much  of  the  numerical  in¬ 
stability  exhibited  by  algorithms  that  are  designed  to 
work  in  an  infinite-dimensional  setting.  This  paper 
attempts  to  generate  some  new  insight  into  both  the  theo¬ 
retical  and  practical  implications  of  Prony's  ideas. 

Young,  T.Y.,  "Representation  and  Analysis  of  Signals,  Part  XVI,  Repre¬ 
sentation  and  Detection  of  Multiple-Epoch  Signals,"  The  Johns  Hopkins 
University,  Baltimore,  Maryland,  AS-405028,  May  1963. 

A  multiple-epoch  signal  consists  of  several  signals  which 
for  one  reason  or  another,  overlap  one  another.  The  purpose 
of  this  report  is  to  design  a  procedure  for  the  detection 
of  the  individual  epochs  of  the  overlapping  signals  and  to 
represent  them  properly.  The  multiple-epoch  signal  may  be 
corrupted  with  random  Gaussian  noise  of  zero  mean.  The 
individual  signals  are  assumed  to  be  representable  by  a 
set  of  exponential  functions  with  acceptable  error,  and  the 
exponents  of  this  set  of  exponentials  are  assumed  known. 

Any  two  adjacent  epochs  are  assumed  to  be  separated  at 
least  T0  seconds  apart. 

For  the  signal  uncorrupted  with  noise  a  criterion  based  on 
the  error  energy  is  described  which  is  of  theoretical  interest. 


This  cricsrion  is  useful  for  the  detection  of  first  epoch 
only.  To  detect  the  second  epoch,  the  first  signal  is  sub¬ 
tracted  from  the  original  signal  by  means  of  the  "comple¬ 
mentary'  operator  concept.  For  a  noisy  signal,  a  likeli¬ 
hood  ratio  criterion  is  proposed.  The  preliminary  experi¬ 
mental  results  carried  out  on  a  digital  computer  justify  the 
theoretical  study. 

Young,  T.Y.  and  W.H.  Huggins,  "Representation  and  Analysis  of  Signals 
Part  XIII:  Discrete  Orthonormal  Exponentials,"  The  Johns  Hopkins  Univer¬ 
sity,  Department  of  Electrical  Engineering,  Baltimore,  Maryland,  AS- 
405028,  May  1963. 

In  dealing  with  the  problem  of  signal  analysis,  one  set  of 
very  useful  component  functions  is  the  orthonormal  expon¬ 
ential  functions.  Previous  studies  have  concentrated 
mostly,  on  continuous  exponentials  with  real  exponents. 

With  recent  progress  in  high-speed  digital  computer,  it 
is  desirable  to  perform  the  signal  analysis  on  digital  com¬ 
puter.  This  requires  the  transformation  of  orthonormal 
exponentials  from  the  continuous  form  into  the  sampled  data 
form,  that  is,  the  z-transformation. 

This  paper  deals  with  digitizing  the  orthonormal  expo¬ 
nentials  with  complex  exponents.  The  z-transforma  of 
continuous  orthonormal  exponentials  are  not  themselves 
orthogonal  in  z-domain.  This  is  due  to  the  fundamental 
property  of  sampling  that  a  signal  with  a  frequency  spectrum 
higher  than  the  sampling  frequency  is  not  exactly 
reproducible  from  its  sampled  data.  For  this  reason,  the 
Kautz  process  is  introduced  in  the  z-domain  to  define  a  sat 
of  discrete  orthonormal  exponentials.  These  discrete 
orthonormal  exponentials  approach  the  ordinary  continuous 
orthonormal  exponentials  in  the  limit  as  the  sampling 
interval  T  approaches  zero. 

Young,  T.Y.  and  W.H.  Huggins,  "On  the  Representation  of  Electrocardio¬ 
grams",  IEEE  Transactions  on  Bio-Medical  Electronics,  Vol.  BME-10,  No.  3, 
pp.  86-95,  July  1963. 

In  dealing  with  the  sophisticated  statistical  analysis  of 
medical  signals  such  as  the  electrocardiogram  (ECG) ,  the 
first  problem  one  encounters  is  how  to  describe  each  ECG 
by  a  few  numbers.  This  is  the  problem  of  efficient  re¬ 
presentation  of  signals,  i.e.,  to  approximate  the  signal 
with  the  smallest  number  of  basis  signals  while  preserving 
the  accuracy  of  the  approximation.  This  paper  begins  with 
a  discussion  of  signal  representation  in  general.  The  con¬ 
cept  of  signal  space  is  introduced,  which  is  very  helpful 
in  understanding  the  ideas  of  signal  representation.  This 
portion  of  the  material  is  of  tutorial  nature. 
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Attention  is  then  directed  to  a  set  of  basis  components 
which  have  found  to  be  very  efficient  for  ECG  representa¬ 
tion,  These  components  are  the  so-called  orthonormal 
exponential  signals.  An  iterative  process  is  developed 
which  enables  one  to  fin-*  a  set  of  mAtched  exponents  for 
the  representation  of  al.  '  >  With  six  pairs  of  such 
exponentials ,  the  average  ►  jf  of  ECG  representation  (QRS 
and  T  waves  only,  leaving  out  P  wave)  is  in  the  vicinity 
of  five  percent.  Experimental  results  will  be  shown. 

Utiing  this  representation,  further  statistical  analysis  may 
be  carried  out  with  ease. 
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